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Abstract. 

The present review is devoted to the muon magnetic moment and its role in 
supersymmetry phenomenology. Analytical results for the leading supersymmetric 
one- and two-loop contributions arc provided, numerical examples are given and the 
dominant tan/3sign(^)/MgygY behaviour is qualitatively explained. The consequences 
of the Brookhaven measurement are discussed. The 2cr deviation from the Standard 
Model prediction implies preferred ranges for supersymmetry parameters, in particular 
upper and lower mass bounds. Correlations with other observables from collider physics 
and cosmology are reviewed. We give, wherever possible, an intuitive understanding 
of each result before providing a detailed discussion. 
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1. Introduction 

The magnetic moment of the muon is one of the most precisely measured and calculated 
quantities in elementary particle physics. It has been measured recently at Brookhaven 
National Laboratory to a precision of 0.54 parts per million (ppm) [1-6]. The Standard 
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Model (SM) theory prediction has reached a comparable level. Thanks to this fantastic 
precision the comparison of theory and experiment is not only a sensitive test of all SM 
interactions but also of possible new physics at the electroweak scale, which is typically 
suppressed by a factor (m^/Mnew ph.)^ and can be expected to contribute at the ppm- 
level. 

Indeed, the experimental result [3] pubhshed in 2001 showed a 4 ppm deviation 
with a statistical significance of almost 3(7 from the SM prediction. This caused a lot 
of enthusiasm and could be nicely explained by a variety of new physics scenarios at 
the weak scale, in particular by weak-scale supersymmetry (SUSY). In the meantime 
several errors in the SM prediction have been corrected, but due to smaller error bars 
on the theoretical and the experimental side, the current 2 ppm deviation still has a 
significance of more than 2a (if the SM prediction is based on e+e^ data, see below). 
Of course, this deviation does not rule out the SM. However, it is intriguing that SUSY, 
which is widely regarded as one of the most compelling ideas for physics beyond the 
SM, would naturally lead to a deviation of the observed magnitude. 

In the present review we describe the rich interplay between SUSY and the magnetic 
moment of the muon. In section 2 we present the analytical results for the SUSY 
contributions at the one- and two-loop level and discuss their main features. In section 
3 the numerical values of the SUSY contributions as functions of the SUSY parameters 
are analyzed. In section 4 we discuss the implications of the experimental result for the 
possible values of SUSY parameters and the correlations with other obscrvablcs. In the 
remainder of this introduction we review the status of the comparison of experiment 
and SM theory and provide some necessary background on SUSY and on the theory of 
magnetic moments. 

1.1. Comparison of experiment and Standard Model theory 

The magnetic moment /I of a particle with mass m and charge e is related to the particle 
spin S by the gyromagnetic ratio g: 

At tree level, QED predicts the exact result g = 2 for elementary spin-|-particles such 
as electron and muon. Quantum effects from QED loop diagrams, from strong or weak 
interactions, or from hypothetical new particles lead to a deviation 

a^l(g-2), (2) 

the so-called anomalous magnetic moment. The theoretical prediction of the anomalous 
magnetic moment of a lepton with mass m is dominated by the one-loop QED 
contribution, the famous Schwinger term a/2n [7], followed by higher-order QED and 
strong interaction effects. Loop contributions from heavy particles with mass M are 
generally suppressed by a factor /M'^ as explained below in section 1.3. Therefore, 
the anomalous magnetic moment of the muon is a factor {jn^jm^'^ ~ 40 000 more 
sensitive to such contributions than the one of the electron. 
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The anomalous magnetic moment of the muon has been determined to an 
unprecedented precision of 0.54 parts per miUion (ppm) at the Brookhaven g — 2 
experiment E821 [1-6]: 

a^^P = 11 659 208.0 (6.3) x 10"^°. (3) 

This is the first magnetic moment measurement that is sensitive to effects from physics 
at the electroweak scale, and it has the potential to constrain and discriminate between 
models of physics beyond the Standard Model (SM). 

Inspired by the success of this experiment, the SM theory prediction of has 
been considerably refined and scrutinized in the last few years (see e.g. [8-12] for recent 
reviews and references). The SM prediction can be decomposed into QED, hadronic 
and weak contributions. The QED contribution is the largest, but it is theoretically 
well under control. The weak contributions are the smallest SM contributions and 
amount to 15.4(0.2) x 10~^°, but they are relevant at the current level of experimental 
sensitivity. Their structure and numerical value are comparable to the ones of potential 
new physics contributions. 

Currently, the hadronic contributions are the main source of the SM theory 
uncertainty. They can be decomposed according to the hadronic subdiagrams into 
vacuum polarization and light-by-light scattering contributions. The hadronic vacuum 
polarization contributions can be inferred from experimental data on the hadronic e'^e" 
annihilation cross section using a dispersion relation. Thus ultimately the precision of 
these contributions is related to the precision of experimental data, which is constantly 
increasing due to ongoing measurements at Novosibirsk, B factories and the (f) factory 
DA$NE. The alternative method of determining the hadronic vacmim polarization from 
data on hadronic r decays suffers from theory uncertainties that arc difficult to assess and 
is not in perfect agreement with the e+e^-based results. The current relative accuracy of 
these contributions is about 1%, corresponding to 7.2 x 10"^*^ error in a^. The hadronic 
light-by-light contributions cannot be related to experimental data and are notoriously 
difficult to evaluate. Early evaluations had a sign error, identified in [13, 14], which led 
to a seemingly large deviation between the experimental result published in 2001 and 
the then current SM prediction [3]. Current estimates vary between 8.6(3.5) x 10""^° [15] 
and 13.6(2.5) x 10"^° [16]. 

For the purpose of the present review we will use the SM theory prediction given 
in [11], based on the most recent e^e^-based evaluations of the hadronic contributions 
[15,17] (see [18] for another recent evaluation, leading to a similar result): 

af^ = ll 659 184.1 (7.2)^""P°1- (3.5)^^^ (Q g^QED/weak ^ ^Q-IO 

= 11659184.1 (8.0) X 10"^°. 
Thus there is a deviation of about 2 ppm between the experimental result and the SM 
theory prediction: 

Aa^(exp - SM) = 23.9 (9.9) x 10"^°. (5) 

In order to appreciate the result further, we put it into historical context and 
compare it to the situation after the first measurement of at CERN in the 1970's, 
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with the result [19] 

^exp,1978 ^ ggg 349 (gS) X IQ-^^. (6) 

This experiment was the first to be sensitive to hadronic contributions to a^, so it was 
interesting to compare its result to the theory prediction without hadronic contributions 
and to the hadronic contributions individually, 

Aa^([exp,1978] - [SM without had]) = 720 (85) x 10"^°, (7) 
aj^<^'^^^« = 667(81) X 10-1°, (8) 

where the central values and errors from Ref. [19] have been used. Without including the 
hadronic contributions there was a gap of more than 8(T between theory and experiment, 
but this gap was beautifully closed by the hadronic contributions. Hence the CERN 
experiment confirmed the existence of hadronic contributions to a^^ and the correctness 
of the SM with a high significance. 

Likewise, the new Brookhaven experiment is the first to be sensitive to weak 
interaction effects, and it is instructive to compare its result to the SM prediction without 
weak contributions and the weak contributions individually, 

Aa^([exp] - [SM without weak]) = 39.3 (9.9) x IQ-^^, (9) 
a^^ = 15.4 (0.2) X 10-1° (10) 

Again, without including the weak contributions there is a gap of about 4(7, which 
estabhshes the existence of contributions beyond the QED and hadronic effects. 
However, in this case the gap (9) is about 2.5 times larger than the SM weak contribution 
(10). 

The deviation (5) or the difference between (9) and (10) has a significance of about 
2.4(j. This is clearly not sufficient to prove the existence of physics beyond the SM. 
The deviation could be due to e.g. statistical fluctuations in the experimental result of 
itself or the experimental data leading to the prediction of aVac.Poi.^ ^j^g imperfect 
understanding of a^^, or a combination of these effects. 

Nevertheless, the result is tantahzing in view of new physics at the weak scale [20]. 
The existence of new physics at the weak scale has been suspected long before the 
measurement of (3), for reasons related e.g. to the naturalness problem of the SM Higgs 
sector, grand unification, or cosmology and dark matter. Generically, new physics at a 
scale M can be expected to contribute at the order mj^/M^ to a^^, up to some numerical 
prefactors. For M ~ M^y such a contribution could easily amount to 2 ppm, and 
this is just the magnitude of the observed deviation (5). In the following we focus on 
supersymmetry as a particularly well-motivated and predictive idea for physics beyond 
the SM. 

1.2. Relevant properties of the MSSM 

The deviation between experimental and theoretical value of could be due to 
contributions from supersymmetry. SUSY at the electroweak scale is one of the most 
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Table 1. The field/particle content of the MSSM. Only 2nd generation (s)leptons and 
3rd generation (s)quarks are hsted explicitly. The mass eigenstates corresponding to 
the electroweak gauge and Higgs bosons and their superpartners are indicated. 



(s)leptons 



(s) quarks 



Higgs gauge bosons 

Higgsinos gauginos 



SM/THDM 



, I^R-, ■ ■ ■ 



SUSY partners L.^ , Ur-, ■ ■ ■ 



, tR, bji, 



tL 



, tR, bu, 



ni,n2 



Hi, H2 



Xl,2,3,4' Xt,2 



compelling ideas of physics beyond the SM (see e.g. [21-25] for reviews). SUSY is the 
unique symmetry that relates fermions and bosons in relativistic quantum field theories. 
It eliminates the quadratic divergences associated with the Higgs boson mass and thus 
stabilizes the weak scale against quantum corrections from ultra-high scales. SUSY at 
the weak scale also automatically leads to gauge coupling unification, and the lightest 
SUSY particle (LSP) can be neutral and stable and constitutes a natural candidate for 
cold dark matter. Moreover, in contrast to many other scenarios for physics beyond 
the SM, the minimal supersymmetric standard model (MSSM) is a weakly coupled, 
renormalizable gauge theory [26], such that quantum effects are computable and well- 
defined, and it has survived many non-trivial electroweak precision tests [27]. 

The MSSM is the appropriate framework for a general discussion of and SUSY. 
The unknown supersymmetry breaking mechanism is parametrized in terms of a set of in 
principle arbitrary soft SUSY breaking parameters. Specific models of supersymmetry 
breaking can be accommodated within the MSSM by suitable restrictions on these 
parameters. 

The MSSM as the minimal supersymmetric extension of the SM contains all SM 
particles and corresponding SUSY partners, see table 1. In addition it also contains 
a second Higgs doublet with associated SUSY partner; hence the MSSM can actually 
be regarded as the SUSY version of the two-Higgs-doublet model (THDM). Two Higgs 
doublets Til, 2 of opposite hypercharge =f1 are required for the cancellation of chiral 
gauge anomalies caused by the corresponding Higgsinos. Thus the field content of the 
MSSM comprises the THDM fields, including five physical Higgs bosons, sec (31)-(33) 
below, scalar SUSY partners of each chiral SM fermion, called sfermions Higgsino 
doublets Hi^2, and U(l), SU(2) and SU(3) gauginos (called bino, winos and gluinos) B, 
W"^'^, g. Right-handed (s)neutrinos as well as non-vanishing neutrino masses are not 
relevant for this review and are ignored. 

Two central MSSM parameters that are of particular importance for are related 
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to the two Higgs doublets. The first of these is the ratio of the two vacuum expectation 
values, 

tan/3 = — . (11) 

SUSY and gauge invariance require that the doublet TYi gives masses to down-type 
fermions, while li.2 gives masses to up-type fermions. As a result, e.g. the top- and 
bottom- Yukawa couplings in the MSSM are enhanced by factors l/sin/3 and 1/ cos/3, 
respectively. In order to avoid non-perturbative values of these Yukawa couplings, tan /3 
is commonly restricted to the range between about 1 and 50. High values tan P = C(50) 
lead to similar top and bottom Yukawa couplings and are therefore favoured by the idea 
of top-bottom Yukawa coupling unification [28]. 

The second important parameter relating the two Higgs doublets is the /^-parameter, 
which appears in the MSSM Lagrangian in the terms 

fjLHiHi - iiFh,'H2 - I^Fh^Hi + h.c. (12) 

The first term describes a Higgsino mass term, while in the other terms Fh^ ^ are 
auxiliary fields whose elimination gives rise to interactions of 7ii^2 with sfermions of 
the opposite type compared to the Yukawa couplings, e.g. to Ti'j'^L^l? ^'^'^ ^2AiA'/i'- 

In addition the MSSM contains a large number of parameters that parametrize 
soft SUSY breaking. Except where explicitly stated we will restrict the number of 
these parameters by neglecting generation mixing in the sfermion sectors. Furthermore, 
we restrict ourselves to the case of it!-parity conservation, since it!-parity violating 
interactions have not much impact on a^. 

The SUSY particles of particular importance to are the smuons, muon-sneutrino, 
gauginos and Higgsinos since they appear in the SUSY one-loop contributions. At 
higher order, also other sectors of the MSSM become relevant, most notably the third 
generation squarks and the Higgs sector. Since in the MSSM all particles of equal 
quantum numbers can mix and these mixings have an important infiuence on the a^- 
prediction, we briefiy discuss the mixing of the individual sectors in the following. 

The sfermions ^ for each flavour can mix, and the mass matrices corresponding 
to the /l, Jr basis read 

^/ - i m,X, Ml^ ) ' ^'^^ 

where 

ML =m} + m\ j + Ml cos2/5(// - Qh%,), (14) 

M\r = m} + m\ j + Ml cos 2/5 Q^s^, (15) 

Xf = Af - n* {cot (5, tan (5) (16) 

with {cot /5, tan /?} for up- and down-type sfermions, respectively, m/, and 
denote the mass, weak isospin and electric charge of the corresponding fermion; 

= sin^ 9]^ = 1 — M^/Ml, where 9^ denotes the weak mixing angle and Mvk,z 
the W and Z boson masses. The quantities Af are soft SUSY breaking parameters 
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for trilinear interactions of sfermions with Higgs bosons of the form /^.-/R-Higgs. The 
remaining entries uil^r of the diagonal elements are governed by the five independent 
soft SUSY-breaking parameters for each generation: 

'mL,i^rn^l = Mq3, m^_^ = M^g, mR^i = Muz, (17) 

mL,ii = mL,i,^ = Mi2, ruR^jx = Mr2. (18) 

We have given these relations for 3rd generation squarks and 2nd generation sleptons 
as these are most important for our purposes. Analogous formulas hold for the other 
generations. The mass matrices can be diagonalized by unitary matrices in the form 

U^'mJU^'^ = diag(m| , m| ) , (19) 

and sfermion mass cigcnstatcs can be defined by 

(I). ./(I). (20, 

The mass of the sneutrino is given by 

ml-ml^ + ^Mlcos2(3 (21) 

and similar for the other generations. 

The superpartners of the charged gauge and Higgs bosons also mix, and the mass 
and mixing terms can be easiest expressed in terms of the Weyl spinor combinations 
-0" = {W~,H{), ij)^ — (W'^^H^)- The mass term for these fields is given by 
il)~Xil)~^ + h.c. with the mass matrix 

where M2 is the soft SUSY breaking parameter corresponding to the SU(2) gaugino 
mass. This matrix can be diagonalized using two unitary matrices C/, V in the form 

U*XV-^^ diag(m^± , ) . (23) 

Similarly, the mass matrix Y corresponding to the superpartners of the neutral gauge 
and Higgs bosons in the basis — [B, W^, Hi^H^) is given by 

/ Ml —MzSw cos 13 Mz-sw sin (3 \ 

M2 Mzcw cos f3 -Mzcw sin f3 

—MzSw cos P MzCw cos P —/jl 

V MzSw sin (3 —MzCw sin (3 —yL J 

where Mi is a U(l) gaugino (bino) soft SUSY breaking parameter and cw — Mw/Mz- 
It can be diagonalized with the help of one unitary matrix N in the form 

N*YN-^ = diag(m^o, . . . , m^o). (25) 

The mass eigenstatcs corresponding to the charged and neutral gauginos and Higgsinos 
are called charginos and ncutralinos, and they are related to the interaction eigenstates 



Y = 



(24) 
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by 

Xt = V.^i^j, (26) 

XT = U,,i,-, (27) 

X- = (28) 

The gluinos do not mix, and their tree-level mass is given by IM3I in terms of the SU(3) 
gaugino mass parameter M3. 

The SUSY parameters /x, Af, Mi_2,3 can be complex. However, not all complex 
phases can appear in observables. The only physical phases of the MSSM (beyond the 
phase in the CKM matrix) are the ones of the combinations 

liAf, /xMi,2,3. (29) 

Hence the frequently adopted convention that M2 is real and positive constitutes no 
restriction. Below, we will adopt this convention only in section 4 and remain general 
elsewhere. 

After spontaneous symmetry breaking the two MSSM Higgs doublets lead to 5 
physical Higgs bosons and 3 unphysical Goldstone bosons. Parametrizing the two 
doublets in the form 

the mass eigenstates are given by 

\ ( cos a sina \ / ci)? \ , , 

(31) 

(32) 

(33) 

where the mixing angle a is related to and the mass of the CP-odd scalar ^4° by 





tan2a = tan2/3 ^2,^2 ^ -- < a < 0. (34) 



The physical Higgs degrees of freedom are the light and heavy CP-even scalars /i'^, 
the CP-odd scalar ^4° and the charged Higgs bosons B.^. For ^ Mz the masses 
of and B.^ are of the order M4. The three unphysical Goldstone bosons are 
eaten to give masses to the Z and bosons. 

For the discussion of two-loop contributions to with Higgs exchange it is 
important to note that the muon receives its mass from the doublet Tii. In the case that 
Ma ^ Mz and tan /? is large the heavy Higgs bosons i?", and are predominantly 
composed of Tii-components. In this case they have larger couplings to the muon than 
the hght Higgs boson 
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1.3. Magnetic moment of the muon and chirality flips 

The magnetic moment is a property of the muon in presence of an electromagnetic field. 
In quantum field theory it is related to the muon-photon vertex function Tf^p^p (p, —p', o) , 
which has the covariant decomposition 

u{p')r^pAp{p, -p', q)u{p) = e u{p') [lpFv{q^) + (p + p')pFM{q^) + ■ ■ ■] u{p) (35) 

for on-shell momenta p, p' and spinors u that satisfy the Dirac equation. The charge 
renormalization condition implies Fy(0) + 2m^FM(0) = 1, and rearranging (35) leads 
to a term 

[1 - 2m^FM(0)]ap,g" (36) 



2m^ 



for — > in the covariant decomposition. This term describes the interaction of the 
muon dipole moment with a magnetic field, and the corresponding gyromagnetic ratio 
is given by gi = 2[1 — 2m^FM(0)], or equivalent ly 

Up = -2m^FM(0). (37) 

In practical calculations it is often useful to extract the form factor Fm before evaluating 
the full vertex function with the help of projection operators [29,30]. 

It is important to understand where the generic behaviour a^'"'^^ oc ^^/Mgygy 
comes from and how it can be modified by additional factors. f The l/Mgugy-behaviour 
for SUSY masses of order MgusY refiects the decoupling properties of SUSY. The m^- 
behaviour of a^, or equivalently the relation Fm oc is related to chiral symmetry, 
and crucial modifications by additional factors are possible. The form factor Fm in (35) 
corresponds to a chirality-fiipping interaction between the left- and right-handed muon. 
If the MSSM were invariant under the discrete chiral transformation 

^U)' ^(/J' (38) 

— IJ'R, — i^R 

of the left-handed doublets and right-handed singlets, Fm and thus would vanish in 
the MSSM. Fm is proportional to because the invariance of the MSSM under (38) 
is broken by the muon mass, or more precisely by all terms in the MSSM Lagrangian 
that are proportional to the muon Yukawa coupling. In each Feynman diagram that 
contributes to a^, the /^-chirality has to be fiipped by one of these terms. The main 
possibilities for the chirality fiip are illustrated in figure 1 and are the following: 

• at a /i-line through a muon mass term, contributing a factor rrifj,, 

• at a Yukawa coupling of Tii to fi^ and fi^ or u^, contributing a factor y^, 

• at a /i-line, corresponding to the transition jiL^jiR, contributing a factor rrinXn ^ 
rrip tan j3 fi for large tan (3 from the smuon mass matrix, 

• at a SUSY Yukawa coupling of a Higgsino to /i and fl or u^, contributing a factor 
t Logarithmic factors like logMgusy/"^^) which can appear at the two- loop level, are disregarded here. 
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• 



OC m. 



fJ'R 



oc tan /? /i 



Pi 



Hr oc rrin tan j3 



fiR oc tan /3 



Figure 1. Possibilities for chirality-flips along the line carrying the /z-lepton number. 



The muon Yukawa coupling is given by 



mf,g2 



- r- , (39) 

where g2 = e/sw, and is thus enhanced by the factor 1/cos/? ~ tan/5 for large 
tan jS compared to its SM value. Hence, while all four possibilities are proportional 
to the muon mass m^, the last three are enhanced by a factor tan/5 for large tan/?. 
SUSY contributions to that make use of these enhanced chirality flips are themselves 
enhanced compared to the generic estimate m^/M|usY- 



2. Analytic results 

In a theory with unbroken SUSY, the gyromagnetic ratio of a charged fermion is exactly 
2 and thus = [31,32], that is the SUSY contribution exactly cancels the SM 
contribution. In the MSSM, SUSY is softly broken, and SUSY contributions to were 
studied already in the early 1980's [33-39], when it was still conceivable that SUSY 
particles were significantly lighter than the W and Z bosons. Although very light SUSY 
particles are now experimentally excluded and SUSY contributions to are suppressed 
as l/M|usY; it is still possible that these contributions are significant. In the following 
we first discuss the generic size of the SUSY one- and two-loop contributions on a more 
intuitive level; then we present the analytic results for the most important contributions. 



2.1. How large can the SUSY contributions be? 

Before presenting the full SUSY one-loop contributions to [33-44], it is instructive to 
discuss the dominant parameter dependence on an intuitive level and to obtain useful 
estimates. As shown in section 1.3 the contributions from SUSY particles of a generic 
mass MsusY are of the order rn'^^/M^^^y, and hence suppressed by a factor / M^^^y 
compared to the SM electroweak contributions. 

However, it has been observed in [37] and further stressed and discussed in [42-44] 
that the SUSY contributions can be significantly enhanced if tan/5 is large. Moreover, 
for large tan (3 the sign of the one-loop contributions is mainly determined by the sign 
of the /i-parameter introduced in (12). We will see here that not only the one-loop but 
also the leading two-loop contributions behave in this way. 
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(NX) 




oc tan/3 /i Ml F(Mi, m^^ ^) 



(N2) 




oc tan/5 /i Ml Ml, m 



(C 



2L J 




OC ml tan /? /i M2 M2, M^J 



2LJ 



tR 



)»- ^ _ ^ >7 



cx ml tan/? ^ {m,X,) F^m^^ ^, M^o] 



Figure 2. Five sample mass-insertion diagrams. Vertices and mass insertions 
are denoted by dots, and the interaction eigenstates corresponding to each hne are 
displayed explicitly. The external photon has to be attached in all possible ways to 
the charged internal lines. The one-loop diagrams (C), (Nl), (N2) have been discussed 
also in [44]. The loop functions F in the results are different in the five cases and 
depend on different masses. 
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It is easiest to understand the leading behaviour with the help of diagrams that 
are written in terms of interaction eigenstates, where the insertions of mass and mixing 
terms and chirality flips are explicitly shown [44] . The five diagrams in figure 2 exemplify 
the main enhancement mechanisms. The basic reason for the tan /3-enhancement is the 
fact that the muon Yukawa coupling in the MSSM is larger by a factor 1/ cos j3 ~ tan {3 
for large tan/3 than its SM counterpart. This Yukawa coupling enters the diagrams 
in figure 2 in the vertices where the muon chirality is flipped, i.e. in the couplings of 
the muon to the Higgsino or Higgs boson in cases (C,N2,C2l,^2l), and in the p,L-fiR 
transition, given by (M|)i2, in case (Nl). 

The second important parameter entering all five diagrams is the //-parameter, 
which governs the mixing between the two Higgs doublets. In all cases, the enhancement 
due to this mixing can be traced back to the fact that 7^2 has the larger vacuum 
expectation value and strongly couples to top quarks, while only Hi couples to muons. 

In diagrams (C,N2,C2l) the //-parameter enters via the Higgsino H1-H2 transitions. 
These transitions enhance the diagrams because the following i72-gaugino transitions 
are by a factor v^z : Vi — tan/3 larger than ^i-gaugino transitions. In diagram (Nl) fi 
enters via the dominant part of the smuon mixing. This mass insertion is obtained from 
the F-term FhJ~L2-i see (12), by replacing H.^ by its large vacuum expectation value. 
Finally, in diagram (t2L) the dominant part of the Higgs-stop coupling originates from 
Fh2^i and thus enables Tii to couple with the top- Yukawa coupling. 

The remaining mass insertions in the diagrams provide additional factors of the 
gaugino mass Mi^2 and stop mixing parameter Xt. They are necessary in order to 
obtain an even number of 7-matrices in the fermion line and in order to connect ti and 
^R, respectively. As an illustration, the relevant factors of diagram (C) are given by 

X22 X12 X22 = A* (5*2^^2) M2 = 92 tan /3 /i M2. (40) 

Combining the enhancement factors of all diagrams leads to the estimates given in figure 
2. They all have a similar form, 

^(C,c,l;Ni,N2) ^ ^ /3 // M2;i F, (41) 

a(f-) ocmJtan/3^ (m,X,)F, (42) 

where M2-i corresponds to (C, C2l) and (Nl, N2), respectively. The loop functions F are 
different in the five cases, depend on the masses appearing in the respective diagrams 
and generally behave as F oc M^gy for large SUSY masses. In these formulas one power 
of is due to the o^-Fm relation (37) and gauge couplings have been suppressed. 

Therefore all leading one- and two-loop contributions are approximately linear in 
tan/3, and their sign is given by the sign of /t, together with the sign of M-1^2 or Xf. 
Generally, all diagrams are suppressed by two powers of the SUSY mass scale. Hence 
the basic behaviour of diagrams (C,C2l,N1,N2) is given by 

2 

^(C,c,l;Ni,N2) ^ ^^^^ sign(//M2;i) (43) 
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if all SUSY masses are set equal to a common scale Msusy- However, it is important to 
keep in mind that the relevant SUSY masses are different in the five diagrams. 

In particular, diagrams (Nl) and {t2h) are special because they increase linearly 
with /i, while the all other diagrams are suppressed for large ^ by their yU-dependent 
loop functions F. Likewise, only the one- loop diagrams are sensitive to the smuon and 
sneutrino masses. If these are large, the one-loop diagrams can be suppressed and the 
two- loop diagrams can become dominant. The chargino-loop diagram (C2l) can be large 
if the chargino and Higgs masses are small. 

The discussion of the stop-loop diagram {i2i) is comphcated by the fact that the 
seemingly linear dependence on [mtXt) is cut off by the requirement that both stop 
mass eigenvalues are positive. This diagram is largest for maximal stop mixing, i.e. if 
{iTLtXt) is large but both eigenvalues are positive and m^-^ ^ m^^i and if mi^ and the 
Higgs boson mass are small. In this case, diagram {t2i) has the behaviour 

4^ ^ ^ ^ tan/3 sign(X,), (44) 

where Msusy denotes here the common mass scale of the appearing Higgs boson and the 
lightest stop. Thus the diagram is linearly enhanced by large /i, and its sign is determined 
by sign(Xt). In the following subsections we will provide the exact analytical formulas 
for all these diagrams and also derive the numerical prefactors in the proportionalities 
(43) and (44). 



2.2. One-loop contributions 

Each diagram that contributes to contains one line carrying the /x-lepton number. 
This fact allows to divide the MSSM one-loop diagrams into two classes: 

(a) SM-like diagrams, where the /x-lepton number is carried only by and/or v^. 

(b) SUSY diagrams, where the yU-lcpton number is carried also by ji and/or z/^. 

The diagrams of the first class involve only SM-particles, and they are essentially 
identical in the SM and the MSSM. The only non-identical diagrams involve two 
couplings of physical SM or MSSM Higgs bosons to the muon line. Owing to the 
additional suppression factor m^/M^ such diagrams are entirely ncghgible both in the 
SM and the MSSM. 

Therefore the SUSY one-loop contribution, i.e. the difference between in the 
MSSM and the SM, is given entirely by the diagrams of the second class. They are 
displayed in figure 3 and involve either a chargino-sneutrino or a neutralino-smuon 
loop. In contrast to the diagrams in figure 2 they are written in terms of interaction 
eigenstates, which is more appropriate for an exact evaluation. The diagrams have 
been evaluated in Refs. [33-36] with various restrictions on the masses and mixings. 
These restrictions have been dropped in Refs. [37-40], and exact results have been 
derived. Later, more comprehensive and general evaluations of these diagrams have 
been presented in the context of particular supersymmetric models [41-43] and the 
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Figure 3. The two SUSY one-loop diagrams, written in terms of mass eigenstates. 
The external photon line has to be attached to the charged internal lines. 



unconstrained MSSM [44] (see also [45, 46] for related results on weak dipole moments 
in the MSSM). We present the general result in the form given in [47]: 

,susY,iL^^xVaf, (45) 

with 



a:. =^T.{-T^i\^L\'+\nL\')F,''{^^rn) + ^MnLnL]F^^ (46) 



16^2^1 12m? ' ™' ' 3m? 

t,rn Mm Mm 

«r = 1^ E {t^(\-'\' + \c^\')F.'i-^) + S^Re[cMi^f (X.)}, (47) 

where i = 1...4 and k = 1,2 denote the neutralino and chargino indices, m = 1,2 
denotes the smuon index, and the couplings are given by 

1 

= V2g^NaU^, + y^N^A. (49) 
4 =-92Vki, (50) 
cf =y^JJk2. (51) 

The kinematic variables are defined as the mass ratios = "^'^o/m'j^^, = m^±/m?^, 
and the loop functions are given by 

Ff (x) = 7^-^-:t^[1 - 6a; + 3x'^ + 2x^ - 6a;2loga;], (52) 



:i-x) 
3 

:i~a;l 



'^{x) = - -^[1 -x2 + 2xlogd, (53) 



2 

Ff (x) = -7 [2 + 3a; - 6x^ + a;^ + 6a; log x] , (54) 

(1 — x)"' 

Ff(x) = ^[-3 + 4x-x2-21ogx], (55) 
(1 — xy 

normalized such that F- (1) = 1. The U(l) and SU(2) gauge couplings are given by 
91,2 = €./{cw,sw}, such that the one-loop contributions are of the order a = e^/(47r). 
A class of large two-loop logarithms can be taken into account by the replacement 
a a(MsusY) (see later for more details). 

For discussing the one-loop contributions a^°' it is noteworthy that the terms linear 
in m^o,± are not enhanced by a factor m^o,± / m^ compared to the other terms. Rather, 
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these terms involve either an exphcit factor of the muon Yukawa couphng 7/^ or of the 
combination t^^it^^2/'"^|„5 which in turn is proportional to (M^)i2 and thus to y^. 
Hence, all terms are of the same basic order ^^/Afgugy, and the terms linear in m^o,± 
are enhanced merely by a factor tan f3 from the muon Yukawa coupling. 

It is instructive to close this subsection by deriving a simple approximation of (46), 
(47) for large tan (3 and the case that all SUSY mass parameters in the smuon, chargino 
and neutralino mass matrices are equal to a common scale Msusy- In this case only 
the terms linear in m^o,± have to be considered, and the loop functions F^{x) can be 
approximated by a Taylor series around x — 1. For example, the factors m^±c;^c^F^(xfc) 
appearing in can be approximated as 



-92yn Yl Uk2Vkim. 




4 mr- 



mr- 



sign(/iM2)Xi2. (56) 



Here terms that are suppressed by 1/ tan/9 or Myy/Msu^Y have been neglected. Note 
that the factors on the right-hand side correspond directly to the mass-insertion diagram 
(C) in figure 2 and the approximation (43). The factors appearing in a^" can be similarly 
approximated. Inserting these approximations, one obtains [44] 



9l 



1927r2 MlusY 



327r2 



m. 



-^SUSY 



sign(/xM2) tan/? 



sign(/iM2) tan/3 



1 + 



1 + 



1 



Ml 



w 



ytan/3' Msusy / 
' 1 Mw \ 
tan/3' Msusy / 



(57) 
(58) 



where real parameters and equal signs of Mi and M2 have been assumed. 



2.3. Two-loop contributions 

It is useful to classify the MSSM two-loop diagrams similar to the one-loop diagrams, 
into 

(a) two-loop corrections to SM one-loop diagrams, where the /x-lepton number is carried 
only by fi and/or z/^,. 

(b) two-loop corrections to SUSY one-loop diagrams, where the /x-lepton number is 
carried also by fi and/or p^. 

The first class contains in particular SM-like diagrams with an insertion of a loop of 
SUSY particles, e.g. of i, b or x"^- Such diagrams are particularly interesting since 
they constitute SUSY two-loop contributions that involve other particles and have a 
completely different parameter dependence than the SUSY one-loop contributions. Most 
importantly, these two-loop contributions can be large even if a^^^^'^^ is suppressed. 
The contributions of this class are exactly known [48,49]. 

The SUSY two-loop contributions of the second class involve the same particles as 
the SUSY one- loop contributions (possibly plus additional ones). Hence they can be 
expected to have a similar parameter dependence as a^^^^'^^. The contributions of this 
class are known in the approximation of leading QED-logarithms [50] . 
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2.3.1. Two-loop corrections to SM one-loop diagrams The MSSM two-loop 
contributions of the first class can be decomposed into a SM- and SUSY-part, 

aSM,2L^^SUSY,2L(a)^ (59) 

where a^'^'^^ denotes the SM two-loop contributions. The genuine SUSY contributions 
of this class can be split into four parts: 

^SUSY,2L(a) ^ ^x,2L ^ ^/,2L ^ ^SUSY,ferm,2L ^ ^SUSY,bos,2L_ (gQ) 

The first two terms correspond to diagrams involving a closed chargino/neutralino or 
sfermion loop, respectively. These diagrams are further categorized according to the 
particles coupling to the muon line, 

^X,2L ^ ^iXVH) ^ ^iXVV) ^ ^iXVC)^ ^ = X. I (61) 

Diagrams where one gauge boson and one physical Higgs boson couple to the muon line 
are denoted as {XV H) with V = j.W.Z and H = h^^H^.A^.H^. Diagrams where 
only gauge bosons or unphysical Goldstonc bosons couple to the muon are denoted as 
{XVV), {XVG).X Sample diagrams are shown in figure 4. 

The remaining two terms in (60) correspond to diagrams involving only SM- or two- 
Higgs-doublet model particles and no SUSY particles. These diagrams are different in 
the MSSM and the SM due to the additional Higgs bosons and the modified Higgs boson 
couphngs. aSUSY,ferm,2L (jeuotcs the difference between the MSSM- and SM-evaluation 
of the diagrams involving a SM fermion (i.e. quark or lepton) loop; likewise, a^^^^'^°®'^^ 
denotes the corresponding difference of the diagrams without fermion loop, the so-called 
bosonic contributions. Sample diagrams are shown in figure 5. 

The SUSY two-loop diagrams can be conveniently evaluated by first applying a 
large mass expansion [51], where the muon mass is treated as small and all other masses 
as large. This results in a separation of scales, and all remaining integrals are of one of 
two types. One type are one-scale two-point integrals with external momentum — 
and all internal masses being either zero or equal to m^. The other type are integrals 
where all internal masses are heavy but the external momentum can be neglected. All 
these integrals and the corresponding prefactors can be evaluated analytically [48, 49] . 
In addition to the genuine two-loop diagrams, one-loop countcrtcrm diagrams have 
to be evaluated. These contain renormalization constants corresponding to charge, 
mass, and tadpole renormalization, which are defined in the on-shell renormalization 
scheme [26,52,53]. 

The diagrams of classes {XV H), X — x, f, f can be calculated in an alternative 
way. In these so-called Barr-Zee diagrams [54], a closed loop generates an effective 
^y-V-H vertex, and this vertex can be evaluated first by a one-loop computation. 
By inserting the result and performing the second loop integral one obtains a simple 
integral representation for the full two-loop diagram. Barr-Zee diagrams were first 
considered because they give rise to important contributions to electric dipole moments 

I Diagrams of the form {XHH), {XHG) etc. in which two Higgs or Goldstonc bosons couple to the 
muon line are suppressed by an additional muon Yukawa coupling and can be neglected. 




Figure 4. Sample two-loop diagrams with closed chargino/neutralino or sfermion 
loop, contributing to a^'^'^ and a^'^'". The diagrams are categorized into classes 
{XVH), (XVG) and (XW), where X = x^'° J ■ V = 1,Z,W^ denotes gauge 
bosons, H = h^,H^,A^,H^ denotes physical Higgs bosons, and G ~ G^^° denotes 
Goldstone bosons. See [48,49] for more details on the possible diagram topologies. 




Figure 5. Sample two-loop diagrams involving only SM- or two-Higgs-doublet model 
particles and either with or without fermion loop. These diagrams are different in the 
MSSM and the SM due to the modified Higgs sector, and this difference constitutes 
the SUSY contributions aSUSY,form,2L ^^^^ ^susY,bos,2L^ 



in extensions of the SM (see e.g. Refs. [54-56]). The contributions from particular Barr- 
Zee diagrams of the classes (fjA^), (fjH), (f'yH), {fW^H"^) to were considered in 
Refs. [57-60], respectively. 

In Refs. [48, 49] the numerical results of all two-loop contributions in (60) were 
compared and analyzed in detail, taking into account that the SUSY parameters are 
constrained by experimental bounds on 6-decays, and other quantities. It turned 
out that the numerical values of the various contributions is very different: 

• The by far largest contributions are the ones from the photonic Barr-Zee diagrams 
(xiH) and (f'yH), where H denotes the neutral physical Higgs bosons h^,H^,A^. 
As explained in section 2.1 they are enhanced by a factor tan/5 and, in the case 
of the sfermion loop diagrams, by the potentially large Higgs-sfermion coupling. 
They can have values up to 

a(x7^^)^ a(/7/^) _ 0(10) X 10"^°. (62) 
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Barr-Zee diagrams with Z or W"^ exchange have a similar parameter dependence 
but are typically smaller by a factor of about 3-5. 

• The diagrams of classes (xW) and {fVV) and the corresponding Goldstone 
diagrams {xVG) and (fVG) involve no enhanced muon-Higgs Yukawa coupling 
and thus no tan/J-enhanccmcnt, and they do not involve any other enhancement 
factors. Their numerical impact is tiny. For SUSY masses larger than 100 GeV 
these contributions are smaller than 0.1 x 10"^''. 

• The genuine SUSY contributions to the SM-like diagrams aSUSY,ferm,2L 
QSUSY,bos,2L (jgpgnd only on ^^^^ and are small. Only for Ma < 200 
GeV they can reach lO"^'^, but for larger Ma they are typically below 0.5 x lO"^'^. 

Hence, for the purpose of the present review, we only present the analytical result for 
the dominant contributions from the photonic Barr-Zee diagrams with physical Higgs 
bosons. They can be written as 

+ E Re[A2AS]/sK+ /M|)], (63) 

S=hO,HO 

<'""'^&;Sp^ E E[ E (A;Q')/RelAj-A|]/,-(myM|)]. (64) 
The Higgs-muon and Higgs-chargino coupling factors are given by 

Xf^HO,AO}^i_Sa C^ A^ (65) 
^{hP,HO,AO} ^ "^^^ [UklVk2{Ca, Sa, -Cp} + t/fc2Vfcl{ - S^, Ca, Sp})- (66) 

In the Higgs-sfermion couphngs we neglect terms that are subleading in tan /3 and that 
give rise to negligible contributions to a^: 

misp 

^1,0,^0^ = ^( - /.*{c„, s J + A,{ - s^, cj) (Ulr Ul (68) 

Xf'"""^ = ^( - /^*{ca, s^} + A,{ - s^, cj) (C/^)* (69) 

The loop integral function fps can be given either as a one- dimensional integral or in 
terms of dilogarithms: 

.idxlog^M 2zr . / l + i/NT 

with y = \/\ — ^z. Note that fps{z) is real and analytic even for 2; > 1/4. The other 
loop functions are related to fps as 

fs{z) = {2z - l)fpsiz) - 2z{2 + logz), (71) 
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Figure 6. Sample two-loop diagrams contributing to ' , i.e. involving a 

SUSY one-loop diagram. The external photon can be attached to all charged internal 
lines, (a) shows a diagram with additional photon loop, giving rise to large QED- 
logarithms. (b) shows a diagram of the class computed in [63]. (c) shows a diagram 
with an additional fermion/sfermion loop. 



f^[z) = '-[2 + \ogz-fps{z)]. (72) 

We remark that useful numerical estimates for the leading two-loop contributions 
can be obtained by taking into account only the tan /3-enhanced terms in the couplings 
and by approximating the loop functions. In the case of the sfermion-loop contributions, 
simple approximations for the most important Barr-Zee diagrams (tHV), (bHV) with 
stop or sbottom loops can be derived [48]. The approximation for [tHV) agrees with 
the estimate (44) discussed in section 2.1, and the approximation for {bHV) has a 
similar form. In the case of the chargino/neutralino-loop diagrams the parameter 
dependence is simpler, and the diagrams with W- and Z-exchange can be included 
in the approximation [49]. We will collect all these approximations below in equations 
(81)-(83). 

2.3.2. Two-loop corrections to SUSY one-loop diagrams The two-loop corrections to 
SUSY one-loop diagrams can be decomposed into two pieces, 

^SUSY,2L(b) ^ ^SUSY,2L(b) AW ^ SUSY,2L(b)_ ^^3^ 

The first piece contains the large logarithm of the ratio Msusy/''^m' "where MgusY is the 
generic SUSY mass scale, and the second piece contains at most small logarithms of 
ratios of different SUSY masses. Sample diagrams are shown in figure 6. The diagrams 
all involve the same particles and the same couplings as the SUSY one-loop diagrams 
(possibly plus additional ones). Hence the overall parameter dependence of o,susY,2L(b) 
and of flSUSY.iL expected to be similar, up to the additional two-loop suppression 

of aSUSY,2L(b) rjnj^g large logarithm is the most relevant enhancement factor. 

The decomposition (73) is analogous to the one of the bosonic two-loop 
contributions in the SM, which have been evaluated in [49,61,62]. In the case of the SM, 
the term enhanced by log M^/m^ is roughly a factor 10 larger than the non- logarithmic 
piece. 

In the case of the MSSM, the logarithmic term is known [50], but for the non- 
logarithmic remainder only a particular subclass of diagrams has been computed recently 
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[63] . In the following we first discuss the logarithmic term, which can be assumed to be 
dominant like in the SM, and then we describe the computation of [63]. 

The large logarithms in (73) are Q ED- logarithms and arise from two-loop diagrams 
that involve a SUSY one-loop diagram and an additional photon loop. The loop integrals 
of such diagrams have an infrared singularity in the limit — > and therefore give 
rise to terms oc logm^. As shown by [50], the appropriate framework to evaluate these 
logarithms efficiently is the framework of effective field theories. 

The relevant effective field theory is obtained from the MSSM by integrating out 
all fields of mass > Msusy and retaining only the muon and photon. All further light or 
heavy SM fields are irrelevant in this analysis and can be ignored. The resulting theory 
is QED with additional higher-dimensional terms, described by 

C,n^-2V2G^Y.Ci{lAOi (74) 

i 

in the notation of [50]. The Oi are higher-dimensional operators. The analysis of [50] 
shows that in the MSSM, like in many new physics models, only one higher-dimensional 
operator has to be considered, namely the one corresponding to the muon anomalous 
magnetic moment, 

H^ = -^m^l^^"'l^Fvp- (75) 

The prcfactors Cj(/x) are renormalization-scale dependent Wilson coefficients, which 
can be determined by matching the effective theory to the full MSSM at the high scale 
Msusy- Determining the Wilson coefficient C//^(Msusy) thus corresponds to the one- 
loop computation of a^^^^. 

By construction, the large logarithms are identical in the full MSSM and the 
effective theory. However, in the effective theory the logarithms can be obtained simply 
from the one- loop renormalization-group running of the Wilson coefficient Ch^ (a*) from 
// = Msusy down io ^ — 171^. This running is described by 

ChM = ChMs^sy) - l{H,, H,) ^ log (Msusy). (76) 

where ^{H^, 11^) is the anomalous dimension of Hfj,. On a diagrammatic level, the 
correspondence of this formula to the two-loop computation in the full MSSM is easy 
to see. In the MSSM, the logarithms arise from diagrams like the one in figure 6 (a). 
Corresponding diagrams in the effective theory arc obtained by contracting the insertion 
of the SUSY one-loop diagram to a point. The resulting diagrams are one-loop 
contributions to if^, involving the effective vertex Their UV- divergence, and thus 
their log/x-terms, determine the anomalous dimension 7(if^,if^). 
The value of the anomalous dimension is 



j{H„H,) = 16. (77) 
As a result, the QED- logarithms in the two- loop contributions to a^^^ are given by 
^susY,2L(b) ^ _4« Msusy ^susy,il ^ susY,2L(b)_ .^g. 
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This logarithmic correction is negative, and it amounts to —7% ... — 9% of the 
SUSY one-loop contributions for Msusy between 100 and 1000 GeV. This result can be 
compared to the case of the bosonic SM two-loop contributions, where the logarithms 
amount to —19% of the SM electrowcak one-loop result. 

As mentioned before, the non-logarithmic terms in a^^^^'^^^'^-' , ^^usY,2L(b)^ 
not known so far. A first evaluation of a subclass of diagrams has been carried out 
in [63]. The considered diagrams involve only sleptons, charginos and neutralinos, and 
only topologies as in figure 6 (b) are taken into account that contain no self-energy 
subdiagrams. These diagrams constitute a finite contribution to (^'^^''^^^^\ However, 
the result of [63] can only be viewed as intermediate because there are more diagrams, 
e.g. containing self-energy subdiagrams or W- or Z-exchange, that would involve exactly 
the same coupling constants, such that non-trivial cancellations could be possible. 

Nevertheless, the investigation of [63] provides a first insight to the possible values 
of the remaining two-loop contributions. The numerical values found in [63] are 
surprisingly large. In a range of SUSY parameters with SUSY masses of the order 
300. . . 500 GeV, the values are mostly below 10~^° but can become up to 2 x 10"^*^, 
which is significantly larger than the corresponding non- logarithmic terms of the bosonic 
SM two-loop contributions. 

2.4- Summary of known contributions and error estimate 

To summarize, the SUSY contributions to up to the two-loop level, i.e. the difference 
of in the MSSM and the SM, are given by 

^SUSY ^ ^SUSY,1L MsusA ^ ^ ^(/,^) 

_^ aix{W,Z}H) _^ ^{f{W,Z}H) _^ ^SUSY,ferm,2L _^ ^SUSY,bos,2L _^ _ _ _ /ygx 

fjj fX fJi fJj ' V / 

where the terms in the first hue have been given analytically in (45), (63), (64), (78). 
Note that the discussion of the two-loop QED-logarithms and equation (76) also show 
that the one-loop result should be parametrized in terms of the running a(AfgusY)- In 
practice, it is sufficiently accurate to approximate a^M^vsv) by a{Mz) = 1/127.9, and 
we define Msusy in the logarithm as the mass of the lightest charged SUSY particle. 

For many applications it should be sufficient to take into account the terms in the 
first line. The explicitly written terms in the second line are known, and they can be 
up to 0{1) X 10~^°, but in the largest part of the MSSM parameter space they are 
much smaller. The dots denote the known but negligible contributions of the type 
(xW), (fVV), the contributions evaluated in [63], and the remaining unknown two- 
loop contributions. 

Handy approximations for the dominant terms are given by 

^susY,iL ^ 13 ^ iQ-10 fl^^^y tan/3 sign{^^M,), (80) 



^SUSY 

al^V") 11v,r,-.o Aang^, ^ 100 GeV ^' 
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._i3xlO-(if?)g)(^)s,.(X,, (82) 

„<M, ._3..,-(!|£)(-i^)(^) (S3, 

The first two are valid if all SUSY masses are approximately equal (note that the 
relevant masses are different in the two cases), and the third and fourth are valid if the 
stop/sbottom mixing is large and the relevant stop/sbottom and Higgs masses are of 
similar size. The result for the iJr)H) contribution has not been here discussed before, 
but it can be understood in the same way as the {t'yH) result [48]. 

In the following we list the missing contributions and estimate the theory error of 
the SUSY prediction (79). 

• Two-loop QED-corrections beyond the leading logarithm (78). The leading-log 
approximation does not exactly fix the scale Msusy in the logarithm and in 
q;(Msusy) (the latter appears in the one- loop result). The exact form of the 
logarithms and of the non-logarithmic terms can only be found by a complete 
computation of the two-loop diagrams with a SUSY one-loop diagram and 
additional photon exchange. The error of the approximation (78) can be estimated 
by varying Msusy in the range 100. . . 1000 GeV to about 2% of the SUSY one- 
loop contributions. If the SUSY contributions to are the origin of the observed 
deviation (5), they are certainly smaller than roughly 50 x 10~^°, and then this 
error is below 1 x 10~^°. 

• Further electroweak and SUSY two-loop corrections to SUSY one-loop diagrams. 
These corrections include two-loop diagrams similar to figure 6 (a) but with 

Higgs- instead of photon-exchange, and like in figure 6 (b) with purely SUSY 
particles in the loops. Given the result for the subclass evaluated in [63], we 
assign an error of ±2 x 10~^° to these diagrams. Note that this is a factor of 10 
larger than the known result of the corresponding non-logarithmic bosonic two-loop 
contributions in the SM [49,61,62]. 

• Two-loop corrections to SUSY one-loop diagrams with fermion/sfermion- loops (see 
figure 6 (c) for an example). This class of diagrams involves in particular top/stop- 
and bottom/ sbottom-loops, which are enhanced by the large 3rd generation Yukawa 
couplings. We estimate the numerical value of these diagrams to be less than 
±0.5 X lO"^*^ for not too hght SUSY masses for the following reasons. SUSY relates 
these diagrams to SM diagrams with top/bottom loops, which amount to about 
0.6 X 10~^° but are not suppressed by possibly heavy SUSY masses. SUSY also 
relates the fermion/sfermion- loop diagrams to pure sfermion-loop diagrams such as 
aj/^^), see (8 2), (83). However, this relation should be most accurate for rather 
small At^b and ji, since the fermion/sfermion-loop diagrams are not enhanced by 
making these parameters large. In that case, the approximations (82), (83) lead to 
values below 0.5 x 10"^'^. 

• Three-loop contributions. In general, three-loop contributions can be expected to 
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be significantly smaller than the two-loop contributions. Two potential exceptions 
are three-loop diagrams that correspond to the two-loop contributions of the types 
(X'yH), {f'~fH) with subloop-corrections to the Higgs-boson masses or the 6-quark 
Yukawa coupling. It is well-known that the one-loop corrections to the Higgs-boson 
masses, in particular to Mh, and to yi, can be very large. Hence, in cases where 
the diagrams with /i-exchange and/or sbottom loop are very large, the missing 
three-loop contributions could amount to x 10~^°. Fortunately, however, the 
influence of the hghtest Higgs boson mass and yb on the (xiH) , (f^H) diagrams 
is small in the largest part of the MSSM parameter space. Hence we neglect the 
theory error associated with the missing three-loop contributions. 

To summarize, we estimate the theory error associated with (79) to 

54^3^ (unknown) = 0.02 af^^^'^^ + 2.5 x 10"^°, (84) 

where the errors associated with the individual classes of missing diagrams have been 
added linearly. If approximated by only the first line of (79), the error increases 

by the neglected contributions in the second line. An upper limit of these can be well 
approximated by [48,49] 

Saf^^^ {2nd fine) = 0.3 (a^f^^^ + a^/^^^) + 0.3 x lO^^". (85) 

It should be noted that the error estimate is deliberately conservative. The later 
numerical analysis, see e.g. table 2, shows that often already the known two-loop 
contributions are much smaller than 10~^°. In these cases, it is reasonable to assume 
that the theory error due to the unknown higher-order corrections is also much smaller 
than the estimate (84). In any case, the theory error of the SUSY contributions is 
smaller than both the current SM theory error and the experimental uncertainty. 

3. Numerical behaviour of the SUSY contributions 

The SUSY contributions to depend on the MSSM parameters in a complicated way. 
In this section we analyze the parameter dependence of the numerical results in detail. 
For each of the one- and two-loop contributions to we will describe which parameters 
are more and which are less influential and explain the dominant behaviour. We discuss 
the parameter choices for which each contribution can become particularly large and the 
numerical values that can be expected for various typical parameter choices. Finally, we 
provide the values of the SUSY contributions to for the benchmark "SPS" reference 
points [64]. 

3.1. One-loop contributions 

The generic behaviour of the SUSY one-loop contributions to is well described by 
(80). The suppression by l/M|ugY, the enhancement oc tan P and the dependence on the 
sign of fi has been explained in section 2.1 using mass insertion diagrams. This generic 
tan P sign(/x) behaviour is illustrated by figure 7, which shows a^^^"^'^^ for various values 
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Figure 7. ay^° ' as a function of a common mass scale Msusy = |a*I = -^2 = 
= iTT-R,!! for various values of tan/3 and /U > (panel (a)), fi <0 (panel (b)). The 
smaller values of tan/3 correspond to smaller |a^^^^'^^|. The preferred la region (5) 
is indicated in grey. 



of tan/3 and both signs of fi. A common mass scale = M2 = tTiL^fx = mB,.ji = Msusy 
has been chosen, except that Afi is fixed by the GUT relation M1/M2 = 5gl/3g2 and 
= 0. Only for small tan/3 there are visible deviations from this simple behaviour due 
to formally 1/ tan /3-suppressed terms. The figure also demonstrates that qiSUSY.il 
be well in agreement with the observed deviation (5), but the small uncertainty in (5) 
results in stringent constraints on SUSY parameters. In the following we consider the 
tan P sign(//) dependence as understood and fix /i to be positive and tan /3 to be large. 

3.1.1. Dependence on mass parameters In general, a^^^^'^^ depends on tan /3 and six 
mass parameters fi, Mi^2, ^R,ii^ ^fi- In the following we will first concentrate on 

the behaviour of a^^^'*''^^ as a function of the four parameters 

M2, mL,p,, niR^p, (86) 

and fix Ml by the GUT relation M1/M2 = 5gl/3g^ ^1/2 and set = 0. 

We visualize the behaviour of a^^^^'^^ for tan/3 = 50 in this four-dimensional 
parameter space with the help of figure 8. In this figure, M2 is chosen to be the smallest 
of the four mass parameters, and /i, mj^.^, m/?^^ are independently varied in the range 
M2 . . . 5M2. The left and right panels correspond to the choices mi fi = M2 (panel (a)) 
and rriL^ix = 5M2 (panel (b)). The difi^erent colours correspond to the values /x = M2 
(dark blue) and n = 5M2 (light yellow). In all regions, the right-handed smuon mass is 
varied between mR^p, = (1 . . . 5)M2. Keeping these parameter ratios fixed, M2 is varied, 
and the resulting a^^^'^^ is plotted as a function of the mass Mlosp of the hghtest 
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Figure 8. o^^^^'^'^ as a function of the mass of tiie ligiitest observable supersymmetric 
particle MlosPj for tan/3 — 50 and the mass parameters varied in the range 
M2...5M2. The blue (yellow) regions correspond to /i = M2 — 5M2) and 
W-R./i = (1 • ■ -5)^/. In panel (a), the left-handed smuon mass is small, rriL.fi = M2, 
in panel (b) m^,.^ = 5A/2- For each case it is indicated which of the mass-insertion 
diagrams (C), (Nl), (N2) are dominant. 



observable SUSY particle [65], defined here as min(m^j, 'rn^±, '"^xl^- 

The two panels show that the /^-dependence is quite different for small and for 
large mL,ji- This intricate interplay between the parameters has been first studied in [44], 
where it has been shown that the mass- insert ion diagrams of figure 2 provide an intuitive 
understanding of the parameter-dependence. In general, the char gino- diagram (C) 
dominates due to its large numerical prefactor in (58). In special parameter regions, the 
diagrams (Nl) with bino-exchange and (N2) with /iR-exchange can become important. 

• Panel (a), dark blue: In the simple case /i = mi^fi = M2, the chargino diagram (C) 
dominates. There is hardly any suppression for heavy mji^f;^, which enters only via 
the smaller neutralino diagrams. Therefore the corresponding dark blue region in 
panel (a) is very narrow. 

• Panel (a), light yellow: If fi is increased, the chargino and most neutralino diagrams 
are suppressed. For large fi the largest contribution can come from the bino-diagram 
(Nl), which is linear in fi. However, (Nl) is only large if both left- and right-handed 
smuons are light. The large spread of the yellow region in panel (a) corresponds 
mainly to this m^j^-dependence of (Nl). The upper border corresponds to small 
rriR^fi (i.e. to yU = 5M2, itil^^ = mn^ji = M2) and leads to values of aSUSY,iL ^^^^ 
even larger than the ones in the blue region. This shows that the largest a^^^'*''^^ 
for a given value of Mlqsp is obtained for 3> m^^ji^^, M2 via diagram (Nl). The 
values at the lower border are suppressed by the large m^jji = 5M2 and are roughly 
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by a factor 3 smaller. 

• Panel (h), light yellow: The situation is reversed in panel (b) for large mL,ji- In 
this case large n does not lead to an enhancement of diagram (Nl) or any other 
diagram, since then all diagrams are suppressed either by the large ^ or the large 
^L,/i- Consequently the result for /i = mL jx = 5M2 is dominated by diagram (C) 
and is almost independent of mp>^jx. Hence the yellow region in panel (b) is narrow. 

• Panel (b), dark blue: For large ttliji but small fi, however, diagram (N2) with fl^- 
exchange becomes important. It is the unique diagram that is not suppressed by 
the large iriL^p,, but it has the opposite sign and a smaller numerical pref actor than 
diagram (C). The lower border of the blue region in panel (b) corresponds to the 
case of small mn^p,, i.e. to niR^j^ — /i — m^p_/5, where the contributions of (C) and 
(N2) almost cancel. For larger ratios than displayed in the plot, m^p^/mn^p^ > 5, 
a^^^^'^'" can even change sign. The upper border of the same region corresponds 
to the case of large rriR^p,, where (C) dominates and (N2) is suppressed. 

Figure 9 shows how the results are modified for different choices of M2 and tan /3. 
Generally, larger values of M2 lead to a suppression of the results without dramatic 
change of the qualitative behaviour, and smaller values of tan/? lead to a suppression 
of the results due to the almost linear tan /^-dependence. Figure 9 (a) shows the same 
as figure 8 (a) except that M2 has been replaced by 5M2, i.e. M2 is the largest of the 
four mass parameters and the other three are varied between M2 /5 . . . M2 . We find that 
the quahtative features are essentially unchanged. The larger value of M2 leads to a 
suppression by a factor of roughly 3, and the largest contributions are now obtained for 
small fj, — M2/5, corresponding to the dark blue region. Panel (b) shows the same again 
but for tan/3 = 5 (and small M2). Here the results are reduced by almost a factor 10 
compared to the case with tan (3 — 50, as expected. 

The two mass parameters we have not explicitly discussed yet arc and Mi. 
Ml has generally less influence than M2 since the Mi-independent chargino diagrams 
dominate in most of the parameter space. However, as exemplifled by the discussion 
of figure 8 (a) there are situations where the bino-exchange diagram (Nl) becomes 
important. If Mi is not tied to M2 by the GUT relation, it is possible to choose 
M2,// » Ml. In this case one can obtain large contributions to a^^^^'^'" from diagram 
(Nl) even if both charginos are very heavy [47]. 

The dependence of a^^^^'-'^'" on Aj^ arises only via the smuon mixing matrix, where 
Afj^ is accompanied by the term yUtan/5, which is typically much larger. Therefore, ^4^ 
is quite insignificant for the SUSY contributions to a^. The small infiuence of was 
verified in [65] by means of a scan of the SUSY parameter space, where 1^4^ | was varied 
up to 100 TeV. 

3.1.2. CP violation and flavour violation So far, we have neglected the possibilities 
of complex phases and generation mixing in the SUSY parameters. The influence of 
complex phases has been studied in [47,66,67]. Clearly, the phase of the //-parameter has 
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Figure 9. a^^^^'^'^ as a function of tiie mass of tiie lightest observable supersymmetric 
particle A^losp, for the same parameters as in figure 8 (a), except M2 5M2 (left 
panel) , tan /3 = 5 (right panel) . 



the most significant impact, corresponding to the sign(yu)-dependence in the real case. 
At the same time, this phase is strongly restricted by negative results for electric dipole 
moment (EDM) measurements. Nevertheless, non-negligible effects of the complex 
phases can be obtained [66,67] even if only small phases that do not violate the EDM- 
bounds are considered. On the other hand, CP invariance is not a critical symmetry 
for the magnetic moment, in contrast to e.g. chiral invariance. Therefore, CP-violating 
phases do not lead to new enhancement factors, and the largest SUSY contributions are 
naturally obtained for real parameters [47]. 

The situation is different for generation mixing in the slepton sector. If the smuons 
can mix with staus, it is possible to obtain contributions where the chirality is flipped at 
a stau-line instead of a smuon-line and thus contributes a factor instead of [44]. 
An example is provided by diagram (Nl), where the jXi^fiR insertion is replaced by a 
tl-tr insertion. Due to the large enhancement factor mr/m^, even a small smuon-stau 
mixing can lead to sizable corrections to qSUSY.il 

3.2. Two-loop contributions 

3.2.1. Chargino /neutralino-loop contributions The two-loop contributions from closed 
chargino/neutralino loops have the same linear tan /?- dependence and 1/M|usy" 
suppression as the one-loop contributions, see (81). The detailed parameter-dependence, 
however, shows interesting differences, as shown in [49]. 

In contrast to the one- loop contributions, the two- loop chargino/neutralino 
diagrams are independent of the smuon mass parameters. They only depend on /i. 
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Figure 10. Panel (a): 
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as a function of the mass of the lightest observable 



supersymmetric particle Mlosp =niin(77i^± , ti^o), for tan/3 — 50 and Ma — 150 GeV 
and for three different ratios /i/M2 = 1, 3, 5 (from top to bottom). Panel (b): Contours 



of 



^SUSY4L 



(xVH),(xVV) 



(solid border) and a' 



SUSY,1L 



alone (dashed line) in the (//- 



M2)-plane for Ma = 200 GeV, mL,^,,^ = 1000 GeV and tan/3 = 50. The contours 
are at (24.5,15.5,6.5,-2.5,-11.5,-20.5) x lO^^*'. Likewise, the colours of the fully 
draws areas correspond to the following values: white: > 24.5 x 10~^°, lightest grey: 
(15.5 . . . 24.5) X lO^i*',. . . , black: < -20.5 x 10"^". The plot has been taken from [49], 
and the contours and regions correspond to Aa^ (exp — SM) as taken in that reference. 



Ml 2, and the Higgs boson mass Ma- Moreover, the dependence on these parameters is 
quite straightforward, while the //-dependence at the one-loop level is very complicated 
due to the different behaviour of diagrams (C), (NX) in figure 2, for example. 

Figure 10 (a) shows the result of the two-loop chargino/neutralino contributions 
a^^^^^'*^^^^^ for three different ratios jji/M2 = 1, 3, 5 as a function of Mlosp, defined 
here as min(m^±, m-^o)- The other parameters are chosen as tan/3 = 50 and M^ = 150 
GeV, and the bino mass parameter Mi is determined by M2 via the GUT relation. For 
j^LOSP > 100 GeV, contributions of up to 5 x 10~^° can be obtained. This is a lot smaller 
than the largest possible one-loop contributions, but it can be a significant correction if 
the one-loop contributions are suppressed by heavy smuons and sneutrinos. 

This can be seen immediately by comparing the two-loop contributions in 
figure 10 (a) with the one- loop contributions in the blue region in figure 8 (b), where 
f^L,ji = 5M2 = 5/i. In this case the two-loop correction can amount to more than 20% 
of the one- loop result. 

The result for a given Mlqsf is largest if /i and M2 are equal. There is no 
enhancement for larger /i similar to the one-loop diagram (Nl). Instead, increasing 
H to n = 5M2 leads to a suppression of the result by a factor 2 ... 3; for even larger 
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ratios than the ones displayed in the figure, the result becomes even smaller. 

For larger values of M^, the result is suppressed as well. We have checked that the 
results in figure 10 (a) would be smaller by a factor 1.5 (for large Mlosp) to 2.5 (for 
small Mlosp) for = 300 GeV instead of 150 GeV. This is in agreement with the 
analysis of the M^i-dependence in [49] for a wider selection of values for fi and M2. 

Figure 10 (a) does not contain the case ji < M2. It is straightforward to show from 
the exact expression (63) that the dominant contribution a^^'^^^ is symmetric under the 
exchange // <-> M2. Therefore the behaviour for // < M2 is very similar to the one for 
/I > M2 and does not need to be analyzed separately. 

The importance of the two-loop corrections is also visible in figure 10 (b), which 
has been taken from [49]. Here the one-loop contributions and the sum of one- and 
two- loop corrections are shown in a contour plot in the (/X-M2) -plane. The smuon mass 
parameters are heavy, mL,R,ii = 1000 GeV, and tan/3 = 50, Ma = 200 GeV. The 
white and lightest grey regions correspond to the Icr-region around the experimentally 
preferred value (the value 24.5(9.0) x 10~^° used for the plot is very similar to the one 
quoted in (5)). For the chosen parameters, this region extends up to /x, M2 < 600 
GeV. However, it is clearly visible that if the two-loop contributions were neglected, 
the contours would shift considerably, and the la region would extend only up to 
M2 < 500 GeV. 

3.2.2. Sfermion-loop contributions The two-loop contributions from sfermion-loop 
diagrams can be even larger than the ones from chargino/neutralino-loop diagrams. 
As discussed in [48,59,60], however, and as the discussion in section 2.1 and formulas 
(82) , (83) show, these large contributions arise only in special corners of parameter 
space. The stop-loop diagrams become large if the stop mixing is large, one stop is 
light, II is very large and Ma is small. Similarly, sbottom-loop diagrams become large 
for large sbottom mixing and if A^/Ma is very large. 

Such parameter regions arc quite restricted for several reasons. Primarily, small stop 
or sbottom masses in conjunction with large tan f3 can have dramatic effects on the SUSY 
predictions for the lightest Higgs boson mass and 6-decays and arc therefore constrained 
by experimental limits on these obscrvablcs. Furthermore, large fi is restricted by 
the requirement that the sbottom mass eigenvalues are positive, and by naturalness 
arguments due to its appearance in the Higgs potential. 

Figure 11 (a) shows the results of the two- loop contributions from sfermion-loop 
diagrams as a function of the lightest sfermion mass in two particular parameter 
scenarios. The parameters have been chosen such that they are in agreement with 
experimental bounds and that the results are sizable. The Higgs boson mass Ma = 300 
GeV, tan/3 = 50, the sfermion mass parameters m^i, mj^j, ttLj^j^, rriL^f, rriR^f have been 
chosen equal to a common scale Msusy- The trilinear parameter At = — 2MgusY; which 
leads to maximal stop mixing and values of Mh, evaluated using FeynHiggs [68-71] 
including higher-order corrections, which arc in agreement with current bounds. 

The light yellow region of figure 11 (a) corresponds to /i — —At — 2MgusY and thus 
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Figure 11. Panel (a): a^J^^ as a function of the mass of the Hghtest observable 
supersymmetric particle Mlosp =niin(TOj^, ra-^^, for the parameter values given in 
the plot. The light yellow (dark blue) region corresponds to [i ~ —At — 2AfsusY 
(/i = — 2Af = 4Msusy), where Msusy is the common sfermion mass parameter, and 
= = 300... 3000 GeV. Panel (b): Maximum values of a^'^^ obtained in a 
parameter scan in the region (87) if the experimental constraints described in the text 
are incrementally applied. The plot has been taken from [48]. 



to a situation where the hghtest stop and sbottom masses are approximately equaL In 
the dark blue region, /i = —2 At = 4Msusy5 and thus the lightest sbottom is significantly 
lighter than the lightest stop. In both regions, Ab = Ar = —300 ... — 3000 GeV. 
The signs of the parameters are such that the sbottom and stop contributions add up 
constructively. The lower borders of the yellow and blue bands correspond to small 
\Ab\ and thus mainly to the pure stop-loop contributions. The width of the bands is 
essentially due to the sbottom-loop contributions, which are approximately linear in Ai,. 

The larger value of /i in the blue region has two effects. On the one hand, the lightest 
sbottom is lighter than the lightest stop, and thus the sbottom-loop contributions have 
more relative importance. Hence the blue band is wider than the yellow band. On the 
other hand, since the stop-loop contributions increase linearly with /x, they are enhanced 
as well in spite of the heavier stops, and hence the blue band lies at higher values of 

The upper borders of the two bands almost exhaust the limits found in [48] for 
the largest possible values of a^'^^. Figure 11 (a), taken from [48], shows these largest 
possible values depending on which experimental constraints are taken into account. 
The allowed parameter range is 

tan/5 = 50, Msusy < 1 TeV, < 3 TeV, 150 GeV < < 1 TeV (87) 

with common sfermion mass parameters m^f = rrij^f: = = rnL,f = Tnn^f = Msusy 
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and = A^. The experimental constraints are > 106.4 GcV, Ap^^^^ < 4 x 10~^, 
BR{B, fj,+ii-) < 1.2 xlO-^ \BR{B X,7)-3.34x 10"^! < 1.5xl0-^ corresponding 
to conservative bounds taking into account experimental and theoretical uncertainties 
(see [48] and references therein). 

If all experimental constraints are ignored, a^'^^ > 15 x 10~^° is possible, in 
agreement with the results of [59,60]. Taking into account the Higgs boson mass hmit 
reduces the maximum contributions drastically, and if all constraints are taken into 
account aj/'^^ turns out to be smaller than 5 x lO"^'^ in the parameter region (87). § 

Significantly larger values of these contributions can be obtained if the sfermion 
mass parameters are non-universal. In particular, if the ratio of ^ and ^ is either 
very large or very small, values of a^'^'" > 15 x 10~^° can be in agreement with all 
experimental constraints [48]. We will come back to this possibility in section 4.6, 
where general scans of the MSSM parameter space are described. 

3. 3. Results for SPS benchmark points 

After having discussed the general parameter dependence of the individual SUSY 
contributions to a^, we present here the results obtained in the Snowmass Points and 
Slopes (SPS) benchmark points [64]. These results provide an overview and reference 
of the SUSY contributions that can be expected in various well-motivated and often 
considered parameter scenarios. Furthermore, we use these results to assess the relative 
importance of the individual one- and two-loop contributions. 

Table 2 shows the results of the known contributions to a^^^^, split up according to 
(79). The QED-improved one-loop and the two- loop contributions with photon exchange 
in the first line of (79) are listed explicitly. The remaining contributions are combined 
into the ones from diagrams with charginos, neutralinos or sfermions and the purely 
SM-like ones. 

The SPS points la, lb, 3, 6 correspond to various minimal supergravity 
(mSUGRA)-scenarios with tan/? between 10 and 30 and SUSY masses in the range 
100 . . . 1000 GeV. They lead to predictions of very close to the observed value (5). The 
same is true for points 7, 8, which correspond to gauge-mediated SUSY breaking. The 
point SPS 2 does not fit so well to (5) because it corresponds to the focus-point region 
in which sfermions are very heavy. SPS 4, 5 involve very large/small tan/?, respectively, 
and therefore yield too high/low values for a^^^^, and SPS 9, which corresponds to 
anomaly-mediated SUSY breaking, involves negative (/xMi 2) and thus leads to negative 

„SUSY 

The two-loop corrections arc very small in all cases. In general, the chargino- 
or sfermion-loop contributions with photon exchange can be the largest two-loop 
contributions, but in all SPS points they are suppressed by moderate values of tan/3 

§ For very low sfermion masses the parameter scenarios used in figure 11 (a) violate the experimental 
bounds, which is why the results displayed in figure 11 (a) are larger than the limits found in figure 
11 (b) at Mlosp ~ 100 GeV. 
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Table 2. Results of the SUSY contributions to in units of 10"^° for the SPS 
benchmark parameter points. The one- loop contributions include the two-loop QED- 
logarithms (78). The SUSY two-loop corrections to SM one-loop diagrams have 

been split up as af,USY,2L(a) ^ ^(xiH) ^ ^(hn) ^ ^SUSY,x+f.re«t + ^SUSY,ferm+bos,2L 

into the photon-loop c;ontributions, the remaining chargino/neutralino and sfermion 
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^SUSY,x-|-f,rest 
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SPS la 


29.29 


0.168 


0.029 


0.056 


0.267 


SPS lb 


31.84 


0.273 


0.044 


0.106 


0.222 


SPS 2 


1.65 


0.032 


-0.002 


0.027 


0.068 


SPS 3 


13.55 


0.078 


0.009 


0.029 


0.187 


SPS 4 


49.04 


0.786 


0.085 


0.288 


0.349 


SPS 5 


8.59 


0.029 


0.135 


-0.046 


0.153 


SPS 6 


16.87 


0.125 


0.015 


0.044 


0.230 


SPS 7 


23.71 


0.236 





0.089 


0.282 


SPS 8 


17.33 


0.163 


-0.001 


0.062 


0.211 


SPS 9 


-8.98 


-0.046 


-0.002 


-0.018 


0.115 



and/or rather high values of Ma- Particularly the sfermion- loop contributions are 
suppressed in addition by the heavy stops and sbottoms in all points except SPS 5. 
The results show, however, that the photon-exchange contributions are larger than 
^susY,x+f,rcst ^ factor ^ 3, in agreement with the discussion in section 2.3.1 and 
the error estimate (85). Likewise, the two-loop corrections from SM-likc diagrams with 
fermionic or bosonic loops are in the expected range ~ (0.1 .. . 0.3) x 10~^°. 

4. Impact on SUSY phenomenology 

In this section we assume that the MSSM is the correct description of physics at the 
weak scale and above and interpret the muon g — 2 measurement within this framework. 
We discuss the impact of the result for on the SUSY parameter space and its relation 
to other relevant observables. 

The deviation Aa^(exp — SM) is positive and larger than the pure SM weak 
contribution, see (5), (9), (10). We have seen before that the MSSM can easily 
accommodate this deviation, preferably for a rather small SUSY mass scale and/or 
large tan (3 and a positive /x-parameter. In the following, we first focus on the general 
MSSM and show that specific, quantitative upper and lower bounds on SUSY masses 
can be derived from a^. Then we compare to 6-dccays, the Higgs boson mass 
and electroweak observables, and ncutralino dark matter. Several of these observables 
exhibit correlations with a^, while others lead to complementary information on SUSY 
parameters. 

Parameter bounds as well as correlations between observables become much more 
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severe in more constrained models than the MSSM. We discuss here the cases of 
minimal supergravity, gauge-mediated SUSY breaking and anomaly-mediated SUSY 
breaking. We conclude the section with a full parameter scan in the general MSSM that 
summarizes the current status of in SUSY. 

Ji^.l. Constraints from on the general MSSM 

Already before the Brookhaven g — 2 experiment, the muon magnetic moment was an 
important observable in SUSY phenomenology. Since the SM theory prediction agreed 
with the older CERN measurement of [19], it was possible to derive lower bounds 
on SUSY masses [33-43]. Even after many of these bounds were superseded by LEP 
bounds, taking into account remained complementary in corners of parameter space 
where light SUSY particles could escape LEP detection [72] . 

In 2001, the publication of the Brookhaven result [3] caused a lot of excitement 
since it showed a deviation of 43(16) x 10~^° (2.6(t) from the SM theory prediction at 
the time (which involved a sign error in the light-by- light contributions). Many authors 
interpreted this result as a possible signal for supersymmetry, and was studied in 
general SUSY models [20,47,65,73-81], with emphasis on correlations and implications 
for particular other observables [82-88], and special scenarios or extended models were 
considered [89-96]. A general conclusion was that significant constraints, in particular 
upper mass limits could be derived even in the general MSSM. 

The current deviation (5) between the corrected SM result and the final 
experimental value is 23.9 (9.9) x 10~^°. This is smaller but statistically almost as 
significant as the one in 2001, and it still allows to stringently constrain the MSSM 
parameter space. Without assuming any specific scenario of SUSY breaking it is possible 
to derive both lower and upperhomids on the masses of SUSY particles. It is particularly 
encouraging that not only one but several SUSY masses can be bounded from above. 

Figure 12 from [97] shows the maximum values of the four lightest SUSY particle 
masses, depending on the value of a^^^^ and tan/3. These values are obtained from 
a scan of the parameters M2, mi,^^, rriR^fi, // up to 2 TeV, assuming the GUT relation 
for Ml. The bounds are independent of the identity of the SUSY particles but they 
mainly restrict the chargino/neutralino and smuon/sneutrino masses as these are most 
relevant for a^^^^. The figure shows meaningful bounds even for the fourth lightest 
SUSY particle. For example, for tan/9 = 10, requiring that a^^^^ > 14 x IQ^^^ 
(a^^^^ > 4.1 X 10^^°), corresponding to a la (2a) band in (5), implies that four SUSY 
particles arc lighter than about 500 GeV (1 TeV) and two of them arc even lighter than 
350 GeV (600 GeV). The upper mass limits are tighter for smaller values of tan/3, but 
all values of tan /3 allow for contributions a^^^^ larger than 20 x 10~^°, so alone does 
not place a lower bound on tan /3 [47, 97] . 

The upper mass limits are very promising for the search for SUSY particles at the 
LHC and ILC. At a linear e'^e" collider, SUSY particles can be pair-produced if they are 
lighter than half of the center-of-mass energy. At the LHC, weakly interacting particles 
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Figure 12. Upper bounds on the masses of the four hghtest SUSY particles as a 
function of fo^ = a^^^^ for tan/3 = 3, 10, 30, 50. The figure has been taken 
from [97], and the dotted vertical lines at a^^^^ = 6, 24 x 10"^" correspond to the 
— Icr lines according to two different SM theory evaluations (using or not using r-decay 
data) available at the time of publication of [97]. Incidentally, they are close to the 
currently preferred value 23.9 x 10^^*^ and the lower 2a contour 4.1 x 10^^". 



like charginos, neutralinos and smuons can be best studied in cascade decays of squarks, 
provided they are lighter than the squarks. In view of figure 12, combined with lower 
squark mass limits from Tevatron, these criteria might be satisfied for a number of SUSY 
particles. 

One might argue that it is too aggressive to require that a^^^"^ is within the la 
or 2a band of (5), especially given the difficulty in assessing the theoretical errors of 
the SM hadronic contributions. Clearly, if e.g. a Sa band is admitted, zero SUSY 
contributions are allowed, and the upper limits on SUSY masses disappear. However, 
interestingly even in a "super-conservative" approach significant parameter constraints 
can be derived [98]. The only requirement made in [98] was that a^^^"^ falls into the 
interval 

-36.8 X 10~^° < af^^ < 89.9 x 10"^° [98], (88) 
which was supposed to be a region that nobody could seriously argue with. At the 
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Figure 13. Lower mass bounds derived from tlie requirement that a^^^^ lies in the 
super-conservative interval (88). For each tan/3 and sign(/i), the region below the 
corresponding line is excluded. The region constrained by direct searches is shaded. 
The figure has been taken from [98] . 



time of publication, it corresponded to the 5cr allowed region, now it corresponds 
roughly to the 6(T region. The bounds derived in [98] from (88) can therefore be 
regarded as definite bounds, independent of any future theoretical or experimental 
developments. In addition to (88), only the following assumptions have been made: 
all SUSY parameters are real, > M2, Mi = M2/2, smuon masses are greater than 
95 GeV, and \ Ay\jmp_^ < 3 in order to avoid electric charge- violating minima. 

Figure 13 shows an example of the (lower) mass bounds that can be derived in 
this super-conservative approach in the plane of the lighter chargino and heavier smuon 
mass. For each tan /3, the region below the corresponding line is excluded. The excluded 
regions are larger for < 0, reflecting the fact that favours positive [i. The excluded 
regions grow with tan /5, but even for small tan /5, excludes regions of parameter space 
that are not excluded by any other experiment. 

General correlations with other observables 

Many observables can be related to in a meaningful way. Here we discuss four 
particularly interesting examples: rare 6-decays, the neutralino dark matter density 
and detection rate, the Higgs boson mass, and electroweak precision observables. Since 
the signs of the parameters play an important role, we fix the convention M2 > for 
simplicity in this and the following subsections. 

We do not discuss in detail the possibility of lepton flavour violation. As mentioned 
in section 3.1.2, /t-f mixing can lead to substantial effects in a^^^^, and conversely the 
measurement of implies bounds on lepton flavour- violating parameters. Furthermore, 
the diagrams contributing to and processes like r ^ fi'y and ^ 67 have a 
similar structure and are correlated. Corresponding detailed studies can be found 
in [91,94,95,99]. 
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4..2.I. B -decays The rare 6-decays b ^ s'y and Bg n'^jj,' are similar to in two 
respects. They are loop-induced, and they involve a chirality flip in the h-s transition. 
Correspondingly, the SUSY contributions to both branching ratios are enhanced by 
tan (3, in the second case even oc tan^ 13. In the case of 6 57 also sign(yu) plays a crucial 
role as it determines whether diagrams with exchange (always positive) and diagrams 
with chargino or gluino exchange (the leading terms are oc sign(yU74j), sign(/xM3), 
respectively) interfere constructively or destructively. The current experimental results 
for the two branching ratios are [100, 101] 

BR(6 ^ S7) = (3.39+°;^?) x 10-^ (89) 
BR(fi, ^ /xV") < 1.0 X 10"^ (95% C.L.). (90) 

Since both results are well in agreement with the SM theory predictions, the possible 
SUSY contributions are bounded from above. As long as minimal flavour violation is 
assumed, this implies for example that destructive interference between the and 
contributions to 6 ^ S7, and thus the specific sign ixAt < is favoured (see [102, 103] for 
reviews and references and [104] for a recent thorough analysis of 6-physics and SUSY). 

The interplay between rare ^-decays and is particularly strong in the framework 
of specific models such as minimal supergravity or gauge-mediated SUSY breaking. Such 
models relate squark and slepton masses and make specific predictions on the signs of 
At and M3 and thus also the sign of /i favoured by 6 ^ 57. Therefore an interesting 
tension can arise between 6 ^ 57 and a^. Both observables favour a particular sign of 
/X, but this might or might not be the same, depending on the values of At and M3. And 
both depend on tan /3 and SUSY masses, but prefers rather large and 6 — > 57 rather 
small SUSY contributions. 

4-2.2. Neutralino dark matter density and detection rate If the lightest neutralino 
is stable it provides a promising candidate for cold dark matter. In this case, two 
observables are of particular interest: the relic density and the detection rate, governed 
by the neutralino-nucleon cross section (for reviews and references see [105, 106]). The 
spin-independent contribution 0-31 to this cross section is yet another chirality-flipping 
quantity, and it is also enhanced by large tan (3. Moreover, for positive this cross 
section is typically larger than lO"^'^ pb, which should be accessible to future detectors, 
while negative /j, would allow for cancellations that could strongly suppress a^i, such 
that neutralino dark matter detection would be out of reach. Hence, the preference of 
for positive /i and not too small tan (3 has significant impact on the prospect of dark 
matter detection [74,78,82,107-110]. Figure 14 from [108] shows the spin-independent 
neutralino-nucleon cross section versus the neutralino mass for a wide range of the 
MSSM parameters /x, M2, tan (3, Ma, A^^t, "m-o, where mo is a common sfermion mass 
parameter. The circles denote points that satisfy the constraint. Imposing this 
constraint increases the minimum value of a^i from 10^^^ pb to 10"^*^ pb, significantly 
improving the prospect for direct detection of galactic neutralinos. 

Under the assumption that standard cosmology (involving a fixed cosmological 
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Figure 14. Neutralino-nucleon cross section versus neutralino mass, for a large range 

of MSSM parameters (sec text), taken from [108]. Only those points are shown for 
which the neutralino relic density could account for all dark matter, but dropping this 
constraint does not affect the conclusions [108]. The circles denote points that satisfy 
the constraint 7 x 10~^° < a^^^^ < 51 x 10~^°, corresponding to a 2a range at the 
time of publication. 

constant and a certain amount of cold dark matter) is valid and that neutralinos 
constitute the only component of cold dark matter, the neutralino relic density VL^ 
is fixed by astrophysical observations, in particular from WMAP [111]. Fixing VL^ 
effectively selects a one-dimensional hyper-surface from the MSSM parameter space. 
However, the dependence of on the MSSM parameters is rather uncorrelated with 
the one of a^. Hence, the measurement of and the increasingly precise determination 
of the cold dark matter density are complementary probes of the MSSM parameter 
space. 

4-2.3. Lightest Higgs boson mass The lightest Higgs boson mass is an important 
observable since it is tightly constrained in the MSSM. At tree level, Mh must be smaller 
than Mz, but the current lower LEP-bound is > 114.4 GeV (95% C.L.) [112]. || Large 
quantum corrections can reconcile the MSSM prediction for Mh with the lower bound, 
but this leads to severe constraints on the MSSM parameter space. The main quantum 
corrections are enhanced by mjlog{mi_^mf^/m^), and large tan/3 and large stop mixing 
can lead to further enhancements (see [27] for a review and references). Therefore, like 

II This bound applies only if the MSSM Higgs boson h is SM-like, which however is the case in the 
largest part of parameter space, in particular for Ma ^ 150 GeV. 
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a^, the M^-constraint favours larger tan/3, but unlike it also favours larger SUSY 
masses, at least in the stop sector. 

Hence there can be a certain tension between the a^- and M/j-constraints, 
particularly in models that connect stop and smuon masses. This tension has become 
stronger recently owing to the latest determination = 171.4(2.1) GeV [113,114]. This 
value is lower than previous ones, and it increases the tendency of the M/j-constraint to 
favour rather large stop masses. In [115], the tension is illustrated for a class of SU(5) 
GUT models by plots of the largest possible a^^^^ as a function of Mh. The higher the 
Higgs boson mass, the lower the possible values for a^^^^. However, the current bounds 
on and can still be simultaneously satisfied in considerable parameter regions. 

4-2. 4- Electroweak precision observables Finally, we briefly comment on electroweak 
precision observables, particularly on M]y and the effective weak mixing angle sin^^eff- 
The present experimental values are [114,116] 

M^P = 80.392(29) GeV, (91) 
sin^^^gP = 0.23153(16). (92) 

The SM predictions for these observables sensitively depend on the values of the SM 
input parameters, in particular on rrit and M^. For the experimentally preferred values 
of the input parameters, the SM predictions agree with quite well and sin^ ^cft- 

It has been noted that the agreement with experiment can be even better in the MSSM 
(see e.g. [117] for an analysis that takes into account the most recent experimental 
data). The SUSY contributions to both observables are dominated by the quantity 
Ap, which is sensitive to the breaking of isospin invariance and thus e.g. to the mass 
splittings in the stop/sbottom sector. These SUSY contributions are enhanced for 
smaller stop/sbottom masses and also depend on the chargino and neutralino masses, 
but they are not particularly sensitive to tan (3. Therefore, as well as the electroweak 
precision observables have a tendency to favour not too heavy SUSY masses, and it is 
fruitful to combine analyses of both kinds of quantities [118-120]. 



4.3. MSUGRA 

In the general MSSM, supersymmetry breaking is parametrized by a large set of free 
parameters. In specific scenarios of supersymmetry breaking these parameters find 
an explanation in terms of some underlying physical mechanism. Typically, then, the 
MSSM parameters can be related to a much smaller set of more fundamental quantities. 
Such supersymmetry breaking scenarios are very predictive, and parameter bounds 
implied by experimental constraints as well as correlations between observables are 
much more stringent than in the general MSSM. 

The first scenario we discuss here is minimal supergravity (mSUGRA) and the 
very similar constrained MSSM (CMSSM) (see e.g. [121] and references therein). In 
this scenario supersymmetry breaking is assumed to take place in a hidden sector 
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Figure 15. (a) The mo-mi/2 plane of the CMSSM for > 0, tan/3 = 10, Aq = 0, 
taken from [132]. The medium pink band and the thin black dashed lines show the 
2cr and la contours for a^. The near- vertical lines correspond to Mh = 114 GeV (red, 
dot-dashed), rn^± = 104 GeV (black, dashed); the medium green area is excluded by 
h — > S7; and the light turquoise band satisfies the relic density constraint. In the dark 
red area the LSP is charged. 

(b) Likelihood function x^^^. for the observables a^, h — *■ 57, Mh, Mw, sin^ ^eff in the 
CMSSM for tan/3 = 10 and various values of Aq. mo is chosen to yield the central 
value of the relic density constraint. This figure has been taken from [120]. 



and to be transmitted to the observable sector via gravitational interactions. The 
Kahler potential of the underlying supergravity theory is assumed to be minimal, i.e. 
in particular to involve no generation- dependent couplings. Furthermore, at the GUT 
scale Mgut ~ 2 x 10^^ GeV the SM gauge interactions unify. The free parameters of 
this model are 

mo,mi/2,v4o,tan/5,sign(/i), (93) 

where mg, "^1/2 and Aq are the universal values of all scalar mass, gaugino mass, and A 
parameters of the MSSM at the GUT scale. The low-energy values of the MSSM soft- 
breaking parameters are determined by renormalization-group running, and the value of 
is determined by the requirement that electroweak symmetry breaking leads to the 
correct value of Mz- (The mSUGRA scenario also leads to a prediction of the gravitino 
mass, but in the present review we will not make use of this prediction.) 

The phenomenology of mSUGRA has been studied extensively, in particular in view 
of the constraint (see [74,77,78,81,88,122-124] for mSUGRA studies focussing on 
and [119,120,125-131] for recent general analyses of mSUGRA). 

As an example, figure 15 (a) from [132] shows the mSUGRA 7710-^1/2 plane for 
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Aq = 0, tan/? = 10, ;U > including contours corresponding to the most important 
obscrvablcs. The region preferred by at the 2a and la level is shown as the medium 
pink band and thin black dashed lines. The Higgs boson mass bound Mh > 114 GeV is 
satisfied to the right of the near- vertical dot-dashed red line at mi/2 ~ 400 GeV, and to 
the right of the black dashed line at mi/2 ~ 150 GeV m^± > 104 GeV. The plot displays 
clearly the tension between M/j, which is increased by larger mi/2, and a^, which prefers 
smaller SUSY masses, but there is a non-vanishing region where both constraints are 
satisfied. For larger tan (5 this overlap region grows. 

The small green region around mi/2 ~ 100 GeV is excluded by the decay h — > 57. 
As mentioned before, the good agreement between SM theory and experiment favours 
destructive interference between the various SUSY contributions. In mSUGRA, the 
dominant contributions are the ones with charged Higgs or chargino exchange, and 
destructive interference requires ^At < 0. Since mSUGRA predicts At < almost 
independently of Aq, positive 11 is preferred by 6 — > 57. It is a non-trivial success 
of mSUGRA that both and b ^ sj prefer the same sign of /x, and this is the 
reason why the green region in figure 15 (a) is so small. However, for larger tan /3 the 
b — > s7-constraint excludes a larger part of the mSUGRA parameter space, and this 
counterbalances the preference of and for larger tan p. Contours for the decay 
Bg — > i.i^/.r are not shown in figure 15, but for large tan/? ^ 50 this decay becomes 
relevant, too. The current GDF limit (90) already begins to constrain the mSUGRA 
parameter space [88,133]. 

The constraint imposed by the requirement that the lightest neutralino relic density 
coincides with the cold dark matter density preferred by WMAP [111] at the 2a level 
is shown by the light turquoise band. This band is very narrow and thus allows to fix 
e.g. mo as a function of mi/2- But since the relic density band is essentially orthogonal 
to the other regions it is possible to satisfy all constraints simultaneously. 

All observables considered in figure 15 (a) have been combined with the electroweak 
precision observables M\y and sin^ ^^cff in [119,120]. In these references, fits within 
mSUGRA have been performed, and remarkably there are mSUGRA parameter points 
that are consistent with all constraints. Figure 15 (b) shows an example of the total as 
a function of mi/2 for tan /? = 10 and various values of Aq (mo is always fixed by the relic 
density constraint). The minimum Xtot — 2-55 is obtained for Aq = mi/2 = 320 GeV. 

In fact, a similarly good fit quality with a Xtot < ^ achieved for all values of 

tan/3 between 10 and 50, and the preferred mass range for mi/2 is always between 300 
and 600 GeV. For lower tan /?, the M^- and a^-constraint are difficult to reconcile. For 
higher tan/?, & ^ 57 is a serious constraint. It should be noted that the significance 
of these tensions has grown recently since the experimental value of the top quark 
mass has gone down [120, 131]. Nevertheless, the fact that mSUGRA fits well to all 
observables and that the preferred mass range mi/2 = 300 . . . 600 GeV is rather low is 
very encouraging also in view of forthcoming SUSY searches at colliders. 
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Figure 16. The Mi -tan /3 plane in GMSB for /i > and for (a) Afmcss = 10^ GeV, 
N^e^s = 1 and (b) M^css = 10^^ GeV, iVmoss = 5, taken from [133]. Tlie siiown 
contours correspond to x 10^" = 10,26,42,58 (blue, dashed), to BR{Bs — > 
(red, solid), and to Mh (black, dot-dashed). The light green region visible only in 
panel (b) is excluded by 6 — > 57. They grey and dark brown regions are excluded by 
mass bounds on SUSY particles and the lightest Higgs boson. 

4.4- Gauge-mediated SUSY breaking 

Gauge-mediated SUSY breaking (GMSB) assumes that supersymmetry breaking is 
mediated from a hidden sector to the observable sector by gauge fields. In the simplest 
case there is an integer number A^'^css of such messenger gauge fields, and these gauge 
fields have mass M^^ss and form vector like 5 + 5 representations of SU(5). A major 
advantage of GMSB over the mSUGRA framework is that fiavour universality naturally 
follows from the symmetries of the messenger interactions (see [134] for a review of 
GMSB). 

The low-energy properties of GMSB are described by the following parameters: 

Mmess, A^mess, A,tan/5,sign(/i), (94) 

where the mass scale A is related to the SUSY breaking scale a/F by A = F/M^ess- A 
determines the overall scale of the soft breaking parameters and can be traded for 
one of them, e.g. Mi. At Mmess boundary conditions are imposed on the MSSM 
soft parameters, and renormalization group running is used to determine the MSSM 
spectrum at the electroweak scale. 

The GMSB predictions for and related observables have been studied in 
[72, 135-137] and compared to the predictions of mSUGRA and other models in 
[47,80,133]. Figure 16 from [133] shows two contour plots in the Mi -tan /3 plane for 
> and for M^ess = 10*^ GeV, N^ess = 1 and M^ess = 10^^ GeV, N^ess = 5. In both 
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cases, the experimental result for can be easily accommodated by GMSB. For any 
given Ml and tan is larger in panel (b) mainly due to the relative suppression of 
the sfermion masses by However, curiously due to the other constraints the 

highest values for a^^^^ can be obtained for smaller A^mcss [47]. 

Like in mSUGRA and b ^ S'j both favour the same, positive sign of /x. Moreover, 
in GMSB both stops and charged Higgs bosons are typically rather heavy, especially 
for small Mmess- As a result, GMSB is not very significantly constrained by the b- 
decays [133] . Conversely, observation of non-standard effects in 6-decays would seriously 
constrain GMSB models. 

4-5. Anomaly-mediated SUSY breaking 

In anomaly-mediated SUSY breaking (AMSB) SUSY breaking takes place in a hidden 
sector and is transmitted to the observable sector via the anomalous breaking of 
superconformal (or super- Weyl) invariance [138. 139]. In this scenario the gaugino and 
scalar mass soft parameters are related to the breaking of scale invariance, expressed 
in terms of the gauge f3 functions and the anomalous dimensions of the matter fields. 
Most notably, the AMSB contributions to the gaugino masses are Mi oc —biOi with the 
one-loop coefficients = (3,-1,-33/5) of the SU(3), SU(2) and U(l) /3 functions. 
While these AMSB contributions to the soft SUSY breaking terms are always present 
in hidden sector models, they are usually subdominant. In the absence of any larger 
contributions, AMSB leads to a very predictive framework that is qualitatively very 
different from mSUGRA or GMSB. 

The parameters of minimal AMSB (mAMSB) are 

Maux,mo,tan/3,sign(/i). (95) 

Maux is a common scale for all soft parameters, and mo is a universal additional scalar 
mass term that does not originate from the super- Weyl anomaly. This term is necessary 
in order to avoid tachyonic sleptons. 

A crucial feature of AMSB is that the signs of M3 and (in the convention 
that M2 > 0) are reversed compared to the mSUGRA and GMSB cases. As a result, 
the 6 — S7 constraint favours negative ^ and thus negative a^^^ . If the observed 
deviation of the experimental from the SM theory value of would have turned out 
to have a negative sign, AMSB would be favoured over virtually all other models of 
SUSY breaking [140,141]. However, since the deviation is positive AMSB is disfavoured 
compared to other scenarios such as mSUGRA or GMSB [65]. This clearly shows how 
combining information on low-energy physics can discriminate between different well- 
motivated SUSY breaking mechanisms. 

Figure 17 shows an example from [133] of the mQ-tan/? plane of mAMSB for 
^aux = 50 TeV and /j, > 0. The region for tan/? > 30 is almost entirely excluded 
by the b ^ constraint (light green). For low tan/3, the Higgs boson mass bound and 
restrict the parameter space. Nevertheless, for moderate values of tan j3 and mo it is 
possible to consistently accommodate 6-decays, and Mh within mAMSB. 
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Figure 17. The mo-tan /3 plane in mAMSB for /i > and il/aux = 50 TeV, taken 
from [133]. The shown contours correspond to x 10^° = 10, 26, 42 (blue, dashed), to 
BR(i?s ji^ (red, solid), and to Mh (black, dot-dashed). The light green region is 
excluded by 6 — s- 57. They grey and dark brown regions are excluded by mass bounds 
on SUSY particles and the lightest Higgs boson. 



4-6. MSSM parameter scan 

In this subsection we present a general, model-independent MSSM parameter scan that 
summarizes the current status of in SUSY. This scan shows the maximum results 
for in the MSSM if all parameters are independently varied and the range of SUSY 
masses for which the MSSM can accommodate the experimental result. The full set 
of one- and two-loop contributions in (79) are taken into account, and the results are 
compared with the current deviation between the experimental and SM result (5). The 
scan can be viewed as an update of the one presented in [65]. 
The MSSM parameters have been varied in the ranges 

|/i|,M2,m^j,m^j, [A/l < 3000 GeV, = 90 . . . 3000 GeV, tan/5 = 50, (96) 

where the upper limit is motivated by naturalness arguments, and the lower limit on Ma 
corresponds to the experimental limit. The parameter tan (3 has not been varied because 
the essentially linear tan j3 dependence of a^^^^ has been sufficiently established before. 
All other parameters have been varied independently, except that Mi has been fixed via 
the GUT relation, = 0, nij^ ^ = rriL^r, = ''^r,t and = A^. These constraints do 
not have much impact on the maximum contributions [65]. The 3rd generation sfermion 
parameters are significantly restricted by experimental constraints on M/j, Ap and h- 
decays. Only parameter points have been considered that are in agreement with these 
constraints, according to the same criteria as in the "weak bounds" of [48]. Note that 
the precise values used in these constraints have not much influence on the maximum 
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Figure 18. Allowed values of a^^^^ as a function of the mass of the lightest observable 
supersymmetric particle Mlosp =niin(TO^± , m^o , ) , from a scan of the MSSM 
parameter space in the range (96) and for tan (3 = 50. The la region corresponding to 
the deviation between experimental and SM values (5) is indicated. The light yellow 
region corresponds to all data points that satisfy the experimental constraints from 
^-decays, Mh and Ap. In the red region, the smuons and sneutrinos are heavier than 1 
TeV. The dashed lines correspond to the contours that arise from ignoring the two-loop 
corrections from chargino/neutralino- and sfermion-loop diagrams. 



results in figure 18 as they affect essentially only the two-loop sfermion contributions. 

Figure 18 shows the results of the scan. The light yellow region corresponds to the 
possible results for a^^^^, given by (79), as a function of Mlosp, if all parameters are 
varied in the range (96). The red region corresponds to the situation that the smuons 
and sneutrinos are heavy, 

m^,_2, m^^ > 1000 GeV, (97) 

while charginos, neutralinos and stops and sbottoms can still be light. A dedicated 
analysis of this situation is of interest since heavy 1st and 2nd generation sfermions 
are sometimes considered as a possible explanation of the absence of observable SUSY 
contributions to flavour-changing neutral currents and CP-violating observables. 

The dashed lines in figure 18 correspond to the results if the genuine two-loop 
diagrams are neglected and only the one-loop contributions and the logarithmic QED- 
corrections (78) are taken into account. 

The yellow region corresponds to all data points and thus to the maximum 
possible values of a^^^"*^ compatible with the parameter range (96) and the experimental 
constraints from 6-decays, and Ap. The SUSY contributions can accommodate the 
observed result (5) within la for LOSP masses below about 600 GeV. For LOSP masses 
below about 440 GeV, the SUSY contributions can be even too large, and thus the a^- 
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measurement significantly restricts the MSSM parameter space in this low-mass region. 

The red region of figure 18 shows that heavy smuons and sneutrinos significantly 
suppress the maximum SUSY contributions to a^. Nevertheless, the contributions can 
be in the la region (5) for Mlosp between about 150 and 470 GeV. 

We close the section with a couple of instructive but more technical remarks. On the 
one hand we analyze which parameter regions give rise to the maximum contributions 
found in figure 18. On the other hand we explain the reasons for the different behaviour 
of the two-loop contributions in the red and yellow regions. 

In the yellow region, corresponding to all data points, the maximum values are 
obtained for /x = 3000 GeV, the upper border of the allowed range, and for M2 ~ mL,R,ji- 
The reason is that the one-loop bino-exchange diagram (Nl) of figure 2 is linear in n and 
therefore dominant for large /i. In comparing our results with the ones of [65], one should 
note that since the maximum contributions arc obtained at the border of the allowed 
multi-dimensional parameter region, an unbiased random generation of parameter points 
as the one used in [65] has difficulties in finding the maximum contributions. Therefore 
the results displayed in figure 1 and equation (3) of [65] slightly underestimate the 
maximum contributions allowed by the employed parameter ranges. Since the maximum 
contributions are obtained for large the two- loop chargino/neutralino contributions 
are negligible, but the two-loop sfermion contributions shift the contour by about 5%. 
Another consequence of the dominance of diagram (Nl) is that the border of the yellow 
contour scales linearly with the maximum value for jj, used in scanning over the SUSY 
parameters. 

The one-loop contributions in the red region are maximized for smuon masses of 
1 TeV and ji — M2. Since therefore /x, M2 and Ma can be simultaneously small, 
substantial two- loop chargino/neutralino contributions are possible. Furthermore, the 
two-loop sfermion contributions can be large as well since the stop and sbottom mass 
parameters are not required to be universal. The largest results in the red region are 
obtained if both rrij^j^ and nij^i are small but is large. In this situation, the 

sbottom-loop contribution can be larger than the results displayed in figure 11 without 
any violation of the experimental constraints on Mh, Ap or 6-decays. 

The significance of the two-loop corrections in the red region is refiectcd in the 
fact that the contour shifts by more than 30% in the low-mass region if the two-loop 
corrections are neglected. 

5. Concluding remarks and outlook 

The final result of the Brookhaven = {g^ — 2)/2 measurement shows a deviation of 
23.9 (9.9) X 10^^°, corresponding to 2 ppm and more than 2o", from the corresponding 
SM prediction based on e'^e~ data. The Brookhaven experiment is the first magnetic 
dipole moment measurement that is sensitive to physics at the electroweak scale, just 
hke the previous CERN experiment was the first to be sensitive to hadronic effects. 
But unlike the CERN experiment, which confirmed the SM prediction of the hadronic 
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effects, tlie Brookliaven experiment would prefer electroweak effects tfiat are about 2.5 
times larger than predicted by the SM. 

The SUSY contributions to are essentially proportional to tan/? sign(/i)/MgusY- 
Hence for moderate or large tan/? these contributions can easily be larger than the 
electroweak SM contributions and thus constitute the origin of the discrepancy between 
the experiment and the SM prediction. Furthermore, in this case strongly favours 
positive which is a very important piece of information for SUSY phenomenology. 
The quahtative behaviour of a^^^^ can be well understood using the mass insertion 
technique. The tan /3-enhancement arises in diagrams where the necessary chirality 
flip occurs at a muon Yukawa coupling, either to a Higgsino or Higgs boson, because 
this coupling is enhanced by 1/ cos (3 ~ tan (3 compared to its SM value. The /i- 
parameter mediates the transition between the two Higgs/Higgsino doublets and 
this transition enhances diagrams because only Hi couples to muons while H2 has the 
larger vacuum expectation value. 

For a quantitative analysis, the MSSM prediction of has to be known with an 
accuracy that matches the one of the SM prediction. This goal has been achieved with 
the calculation of the full one- loop and leading two- loop contributions, as given in (79). 
Since the remaining theory uncertainty is dominated by one particular class of unknown 
two-loop diagrams, an even more satisfactory MSSM prediction could be obtained by 
evaluating this missing class of diagrams. 

If supersymmetry is assumed to exist, many non-trivial details about the SUSY 
spectrum can be derived from a^. For example, for tan/? = 10 four SUSY particles must 
be lighter than 500(1000) GeV at the la{2a) level. And even in a very conservative 
interpretation, where a 5(7 deviation in is tolerated, significant lower bounds on 
the SUSY masses can be derived in a model-independent way. These bounds cannot 
be obtained from any other observable and establish as one of the most important 
indirect probes of SUSY. 

There arc several other observables that are relevant for SUSY phenomenology and 
provide constraints on SUSY parameters, although is unique in its largest discrepancy 
between experimental and SM values. On a qualitative level, the relation between 
and these other observables is the following, favours rather large tan/? and/or small 
SUSY masses and positive /x. If neutralinos constitute a component of dark matter, the 
dark matter detection rate is enhanced by large tan (3 and positive 11. Therefore, the 
result is encouraging for future dark matter detection experiments. 

The preference of for rather light SUSY masses is supported by the experimental 
results for the electroweak precision observables Mw and sin^ ^eff, which favour small but 
non-zero SUSY contributions. The preference of for not too small tan /? is supported 
by the constraint derived from the lower bound on M^. However, the bound on also 
favours large SUSY masses (particularly stop masses) over light ones. Furthermore, the 
MSSM parameter space for large tan /? is significantly constrained by the rare 6-decays 
6 — > 57 and Bs —>■ A*"*"/^"- The decay b ^ shows a similar dependence on sign(/x) as 
a^. This already constitutes a crucial test of minimal models. For instance, in mSUGRA 
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or GMSB, but not in AMSB, and h ^ favour the same sign of /i. 

Given the tension between all these observables it is non-trivial that all constraints 
can be simultaneously satisfied, even in the simplest but well-motivated models such as 
mSUGRA or GMSB. Mainly driven by a^, the parameter points satisfying all constraints 
prefer a number of (though not necessarily all) SUSY particles to be light. This result 
is clearly encouraging for the SUSY search at the LHG and the ILG. 

We have summarized the current status of and SUSY in figure 18, which shows 
the possible values of a^^^^, compared with the observed deviation between experiment 
and the SM prediction. The MSSM parameters have been scanned over in the range 
allowed by all relevant constraints from other collider experiments, and in the evaluation 
of a^^^^ all known one- and two- loop contributions have been taken into account. The 
result confirms again that SUSY can accommodate the result consistently with all 
other constraints and that the preferred mass range is rather low. On a more technical 
level, the figure also shows that two-loop SUSY effects can be important. 

In spite of the impressive current status, progress on is important and will come 
both from the experimental and the theoretical side. Already a small improvement of the 
precision could be sufficient to increase the discrepancy between the experimental and 
SM values of to 4-5(7 and thus to estabhsh the existence of non-SM contributions. 
Gurrently the SM theory prediction has a larger uncertainty than the experimental 
value, and within the SM prediction the main uncertainty is related to the hadronic 
vacuum polarization. However, this uncertainty can be expected to be significantly 
reduced as a result of ongoing measurements of e+e^ hadrons by CMD2 and SND in 
Novosibirsk and, using radiative return, by B factories and KLOE. Another important 
source of theoretical uncertainty is related to the hadronic hght-by-light scattering 
contribution. It is a very challenging theoretical task to evaluate this contribution, 
and current error estimates vary between 2.5 ... 4 x 10~^°. However, further progress in 
the understanding of this contribution can be expected, and one can hope that the full 
SM theory uncertainty can be reduced to below 4 x 10~^° in the foreseeable future [17]. 
Finally, progress can be expected from the experimental determination of a^. Since 
the current measurement precision is statistics limited, a new experiment using similar 
methods with only straightforward improvements could reduce the uncertainty by a 
factor of 2.5 or more, and an according plan has already been outlined [6, 142]. 

After more than 40 years of experimental and theoretical progress, the observable 
a^ii has become a sensitive probe of physics at the electroweak scale. Already today it 
is one of the most important constraints of physics beyond the SM. In the near future, 
particle physics will enter a new era where the detailed structure of physics at the 
electroweak scale and above will be unravelled by experiments at the LHG and possibly 
later the ILG. Not only new particles like smuons or charginos might be discovered, but 
also the other indirect constraints from the observables mentioned above will become 
more stringent. The magnetic moment of the muon will provide a crucial cross-check 
and complement of the forthcoming experiments. 
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